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CHAPTER  I :  INTRODUCTION 


Statement  of  the  Problem 

Extreme  wave  conditions  have  been  the  cause  of  property  loss,  suffering, 
injuries,  and  death  since  man  first  approached  the  sea.  Coastal  engineers 
therefore  have  always  attempted  to  build  works  that  would  withstand,  with  lit¬ 
tle  or  no  damage,  the  worst  impact  of  waves  from  very  rare  events.  The  direct 
effect  of  sea  waves  striking  coastal  structures  has  long  been  recognized  as  a 
critical  phenomenon  with  respect  to  structural  integrity  during  a  storm  at 
sea.  The  hydraulic  impact  of  individual  waves  has  traditionally  been  the  spe¬ 
cific  force  used  as  the  basis  of  structural  design  criteria;  therefore,  char¬ 
acteristics  of  the  worst  few  waves  of  a  hypothetical  extreme  event  have  been 
estimated  for  application  in  most  design  computations.  Rubble-mound  struc¬ 
tures,  constructed  of  layered  quarrystone  or  concrete  shapes  and  built  for 
centuries  as  wave  barriers  (breakwaters  and  jetties)  or  shore  protection  (re¬ 
vetments),  are  usually  designed  in  this  fashion. 

The  limits  of  functional  performance  of  coastal  structures  have  recently 
become  more  critical  with  respect  to  overall  economic  optimization.  Public 
financing  of  coastal  works  has  been  more  difficult  to  arrange  than  in  past 
decades.  The  concept  of  designing  a  structure  to  be  stable  during  a  very  ex¬ 
treme  storm,  but  to  be  less  than  1UQ  percent  effective  in  some  extreme  events 
of  lesser  intensity,  has  been  in  the  minds  of  coastal  engineers  in  an  effort 
to  conceive  affordable  harbor  or  shore  protection  plans.  Life  cycle  cost  also 
is  receiving  much  more  scrutiny,  particularly  with  respect  to  expensive  mobi¬ 
lization  and  challenging  construction  techniques  required  for  repairs  at  many 
coastal  projects.  The  bulwarks  of  extreme  conservatism  in  coastal  engineering 
design  practice  are  beginning  to  buckle  under  pressure  for  more  precise  esti¬ 
mates  of  structural  integrity  and  functional  performance.  These  estimates  may 
someday  approach  the  precision  of  those  now  required  for  design  of  buildings 
and  bridges. 

One  critical  question  in  many  new  optimized  designs  is  "What  is  the  ef¬ 
fect  of  duration  of  exposure?"  Sandy  beaches  commonly  change  their  shapes  to 
a  more  stable  configuration,  given  sufficient  exposure  to  severe  wave  condi¬ 
tions,  in  theory  approaching  a  new  equilibrium  (Bruun  1954).  Some  radical  new 
rubble-mound  concepts  attempt  to  emulate  this  effect  (Delft  Hydraulics 


Laboratory  1985) .  Laboratory  experiments  which  simulate  natural  irregular 
waves  also  have  shown  some  duration  effects  on  rubble  mounds  of  more  tradi¬ 
tional  design  (Graveson  et  al.  1980;  Van  der  Meer  and  Pilarczyck  1984;  and 
Tenaud  et  al.  1981).  The  open  literature  contains  little  specific  guidance, 
however,  for  researchers  or  designers  to  estimate  the  duration  of  a  given 
intensity  of  extreme  wave  conditions. 

Purpose  and  Objectives 


The  purpose  of  this  work  is  to  investigate  the  duration  of  extreme  wave 
conditions  estimated  from  hlndcast  wave  data,  with  a  view  toward  developing  a 
means  to  characterize  the  variation  of  these  durations  for  use  in  design  of 
coastal  structures.  Hlndcast  wave  data,  which  are  discussed  later  in  more 
detail,  are  one  of  the  most  valuable  tools  of  coastal  engineers,  primarily 
because  weather  data  on  which  they  are  based  typically  exist  for  much  longer 
periods  of  record  than  other  wave  information  sources.  The  20-year  (1956- 
1975)  Wave  Information  Studies  (WIS)  database  of  hlndcast  wave  data  prepared 
and  maintained  by  the  US  Army  Engineer  Waterways  Experiment  Station  (WES) 
(Brooks  and  Corson  1984)  is  a  key  source  of  wave  information  in  many  US  Army 
Corps  of  Engineers  projects  since  it  now  extends  along  most  of  the  coastline 
of  the  United  States. 

The  specific  objectives  of  this  study  were  to  (1)  review  existing  liter¬ 
ature  regarding  the  duration  of  extreme  wave  conditions  and  related  topics; 

(2)  formulate  a  practical  means  of  identifying  individual  events  of  extreme 
wave  conditions,  relying  on  the  intensity  of  wave  conditions  as  represented  in 
the  WIS  database  and  associated  publications;  (3)  address  the  probability  dis¬ 
tribution  of  extreme  event  durations  by  fitting  selected  distribution  func¬ 
tions  to  representative  data;  and  (4)  address  the  possible  relation  of  an  ex¬ 
treme  event '8  duration  to  the  peak  conditions  during  the  extreme  event  by 
regression  analysis. 

Organization 

This  report  presents  reviews  of  pertinent  statistical  concepts  and  tech¬ 
niques,  considerations  regarding  the  characterization  of  wave  conditions,  and 
the  specific  nature  of  WIS  hlndcast  data  before  proceeding  to  describe  the 


progress  toward  and  conclusion  of  the  four  objectives  stated  above.  An  over¬ 
all  summary  and  statement  of  conclusions  then  is  followed  by  Appendix  A  con¬ 
taining  figures  and  tables  which  were  not  presented  in  the  main  text  for  the 
sake  of  continuity  and  space  conservation.  Appendix  B  includes  pertinent  wave 
information  transcribed  from  the  WIS  database.  Appendix  C  includes  a  listing 
of  the  computer  program  STRMDIST  which  was  used  to  identify  extreme  events, 
define  durations,  and  fit  parameterized  distribution  functions  to  both  the 
durations  and  peak  wave  heights  of  extreme  events  identified.  Appendix  D  in¬ 
cludes  the  command  file  for  the  commercial  statistical  software  package  SPSS 
(Nie  et  al .  1975),  which  was  applied  to  address  the  relationship  of  extreme 
event  duration  to  peak  wave  conditions. 
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CHAPTER  II:  REVIEW  OF  PERTINENT  STATISTICAL  CONCEPTS 


Continuous  Frequency  Distributions 

The  primary  tools  of  this  study  are  statistical  procedures  which  address 
the  variability  of  parameters  of  interest,  specifically  duration  of  extreme 
events  at  sea  and  their  peak  intensity.  A  brief  review  of  pertinent  statisti¬ 
cal  concepts,  which  are  critical  to  understanding  the  methods  and  conclusions 
of  the  analysis,  is  presented  below. 

Continuous  random  variables  are  variables  whose  values  are  measured  on  a 
continuous  scale,  as  opposed  to  their  discrete  counterparts  such  as  rolling 
dice  or  coin  flipping.  Most  natural  phenomena  of  varying  intensity  as  mea¬ 
sured  by  instruments  are  treated  as  continuous  random  variables.  The  proba¬ 
bility  that  the  value  of  a  particular  random  variable,  x  ,  will  fall  within  a 
certain  range  can  be  estimated  by  application  of  its  probability  density  func¬ 
tion,  f(x)  ,  which  is  analogous  to  a  histogram  for  discrete  variables.  The 
following  two  conditions  apply  in  defining  probability  density  functions: 

f(x)  ^  0  for  all  x  within  the  domain  of  f 

and 

00 

/  f ( x)  dx  =  1  (1) 

_00 

The  probability  that  x  will  fall  within  the  range  from  a  to  b  is  given 
by: 


P(a  *  x  *  b) 


dx 


(2) 


Technically  the  probability  of  x  taking  on  a  value  of  exactly  a  or 
b  is  zero,  but  since  physical  measurements  cannot  be  infinitely  accurate,  the 
interval  from  a  to  b  can  be  considered  inclusive.  A  transformation  of  the 


probability  density  function  into  its  corresponding  distribution  function, 
F(x)  ,  allows  more  expedient  computation  of  probabilities: 


F(x) 


fCt)  dt 


(3) 


where  f(t)  is  the  probability  density  function  of  a  dummy  variable  t  . 

The  value  of  F(x)  varies  between  0  and  1.  The  probability  that  x 
will  have  a  value  equal  to  or  less  than  a  is  F(a)  .  The  probability  that 
x  will  have  a  value  between  a  and  b  is  F(b)  -  F(a)  .  The  corresponding 
probability  density  function  is: 


f(x)  -  (4) 

It  is  important  to  define  the  domain  of  f  and  that  this  domain  include 
all  the  values  of  x  of  interest.  Furthermore,  the  function  f  must  be  in- 
tegrable  within  this  domain  (and  F  differentiable)  for  the  above  definitions 
to  appJy  (Miller  and  Freund  1985). 

Distribution  Parameters 


The  mean  or  expected  value  of  x  is  defined  by: 


U 


xf(x)  dx 


(5) 


The  variance  of  probability  density  function  is  the  expected  value  of 
the  squared  deviation  from  the  mean,  given  by: 


o 


2 


(x  -  v)?  f(x) 


dx 


(6) 


5 


The  variance,  o  ,  and  Its  square  root,  the  standard  deviation,  a  , 
are  both  measures  of  the  spread  of  the  probability  density  about  the  mean. 

The  standard  deviation  is  expressed  in  the  same  units  as  x  and  y  .  A  small 
variance  or  standard  deviation  implies  a  strong  central  tendency  while  large 
values  imply  significant  spread  or  "variance"  of  x  values  (Miller  and  Freund 
1985). 


The  Poisson  Distribution 


A  wide  variety  of  distribution  functions  have  been  formulated  by  re¬ 


searchers  and  statisticians  which  have  been  shown  to  describe  well  the  behav¬ 


ior  of  certain  random  variables  which  occur  in  nature.  One  such  function  is 


the  Poisson  distribution,  defined  by: 


,x  -X 

X  e 


for  x  -  0,  1,  2,, 


This  is  a  discrete  distribution  which  has  important  associations  with  the  con¬ 
tinuous  distributions  that  have  been  applied  to  describe  weather-related  vari¬ 
ables.  Specifically,  the  Poisson  distribution  has  been  applied  to  describe 
the  number  of  occurrences  of  events  taking  place  randomly  over  continuous 
interval**  of  time.  The  parameter  X  is  both  the  mean  and  the  variance  of  the 
Poisson  distribution.  A  key  assumption  behind  application  of  this  distribu¬ 
tion  is  that  the  probability  of  an  occurrence  for  the  type  of  event  in  ques¬ 
tion  during  a  small  interval  of  time  must  not  depend  on  what  happened  prior  to 
that  time.  A  random  process  which  fits  this  criterion  is  called  a  Poisson 


process. 


The  Exponential  Distribution 


A  continuous  distribution  which  is  often  associated  with  the  Poisson 


distribution  is  the  exponential  distribution,  given  by: 


f(x)  -  - 


for  x  >  0  and  8  >  0 


■  0  elsewhere 


The  corresponding  distribution  function  is: 

F(x)  -  1  -  e'x/e  (9) 

The  mean  and  standard  deviation  of  a  variable  represented  by  an  exponential 

2 

distribution  are  both  g  and  the  variance  is  g  .  This  distribution  is 
often  used  with  Poisson  processes  to  model  the  waiting  time  between  successive 
occurrences.  If  the  x  parameter  of  a  Poisson  distribution  is  the  average 
number  of  occurrences  in  time  T  ,  then  the  average  rate  of  occurrences  per 
unit  time  is  \/T  .  The  corresponding  exponential  distribution  parameter  is 
g  *  T/x  .  This  relation  and  the  fact  that  both  distributions  are  fully  de¬ 
scribed  by  a  single  parameter  make  them  easy  to  use  in  a  wide  range  of  appli¬ 
cations  dealing  with  the  frequency  of  and  waiting  time  between  discrete 
events. 

The  Weibull  Distribution 

Another  distribution,  which  is  widely  used  to  model  the  variation  in  in¬ 
tensities  of  natural  extremes  such  as  flood  elevations  and  storm  intensities, 
is  the  Weibull  distribution,  where: 

for  x  >  0  ,  a  >  0  ,  g  >  0 

■  0  elsewhere  (10) 

The  corresponding  Weibull  distrubution  function  is  very  similar  to  the  expo¬ 
nential  distribution: 

F(x)  -  1  -  exp  £-  ^ 

The  parameter  a  is  the  "shape  parameter"  which  defines  the  basic  shape 
of  the  function.  The  g  parameter  is  the  "scale  parameter"  which  determines 
the  degree  of  spread  along  the  abscissa  (Isaacson  and  MacKensie  1981).  The 
mean  and  variance  of  the  Weibull  distribution  are: 
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"  -  er(i  -  i) 


(12) 


(13) 


The  gamma  function  is  given  by: 


r(z) 


(14) 


The  Welbull  distribution  has  two  parameters  which  make  it  actually  a  family  of 
functions.  A  three-parameter  form  Is  sometimes  used  to  provide  further  flexi¬ 
bility  In  adapting  the  distribution  to  certain  phenomena,  where: 


F(x) 


exp 


for  £  >  0 


(15) 


The  parameter  £  Is  a  "location  parameter"  which  locates  the  position  of  the 
probability  along  the  abscissa  (x-axls).  In  the  particular  case  of  the 
Welbull  distribution,  £  Is  in  effect  a  lower  limit  to  values  of  x  .  The  e 
parameter  Is  often  taken  as  zero  in  practice.  The  Welbull  distribution  re¬ 
duces  to  the  exponential  distribution  when  a  -  1  and  £  «  0  (Isaacson  and 
MacKensle  1981). 

The  Rayleigh  Distribution 

The  Welbull  distribution  reduces  to  a  Rayleigh  distribution  when  ®  ■  2 
and  £  ■  0  ,  a  function  widely  used  to  model  the  distribution  of  wave  heights 
passing  a  point  during  a  stationary  sea  state.  The  term  "stationary"  refers 
to  the  common  assumption  that,  for  practical  purposes,  statistical  properties 
of  ocean  waves  tend  to  be  time  Invariant  during  a  period  of  a  few  minutes  to 
an  hour  or  more.  The  time  for  significant  changes  to  occur  in  a  sea  state  Is 
thus  assumed  to  be  substantially  longer  than  the  time  necessary  to  measure  the 
form  of  a  few  hundred  waves  passing  a  fixed  point.  The  Rayleigh  distribution, 
for  this  purpose.  Is  often  expressed  in  the  form: 
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-2(H/H  )2 

F(H)  -  1  -  e  8  (16) 

where  H  Is  an  Individual  wave  height  In  a  sea  state  and  Is  the  "signif¬ 

icant  wave  height,"  also  defined  as  the  average  of  the  highest  1/3  waves. 

This  relation  has  been  found  to  be  quite  accurate  In  most  conditions  at  sea, 
with  the  exception  of  waves  nearing  the  point  of  breaking  In  shallow  water 
(Massle  1976).  The  corresponding  probability  density  function,  mean,  and 
variance  of  this  form  of  the  Rayleigh  distribution  are: 


f  (x) 


u 


0.627  H 

8 


(17) 


(18) 


2 


a 


(a  -  0.779  H  ) 

s 


The  Extremal  Type  I  Distribution 


(19) 


This  distribution',  sometimes  called  the  "Gumbel"  or  "Fisher-Tippet  Type 
I"  distribution,  also  Is  frequently  applied  to  model  natural  extremes  such  as 
storm  intensities  (Gumbel  1958) .  The  probability  density  and  distribution 
functions  have  the  following  forms: 


_e-[(x-e)/*]  _[(x_c)/+] 

f(x)  -  £ - _£ - 


for 


— oo  <  X  <  00 

—00  <  £  <  OO 
6  >  0 


-[  (x-e)/B] 

F(x)  -  e"e 


(20) 


(21) 


The  mean  and  variance  are: 


V  -  e  -  y6 


(22) 


(23) 


where  y  m  Euler's  constant  “  0.5772  .  The  Extremal  Type  I  distribution  is 
also  a  two-parameter  family  of  functions,  in  this  case  with  a  shape  parameter 
of  a  -  1  In  keeping  with  the  usual  practice  for  application  to  weather- 
related  phenomena  (Isaacson  and  MacKensle  1981  and  Andrew  et  al.  1985).  The 
e  parameter  is  again  the  location  parameter  and  6  the  scale  parameter.  The 
Extremal  Type  I  distribution  is  not  constrained  to  positive  values  of  x  . 

Figure  1  illustrates  the  relative  form  of  the  Exponential,  Welbull, 
Rayleigh,  and  Extremal  Type  I  distributions.  The  Exponential  and  Rayleigh 


Cumulative  Probability 

Figure  1.  Relative  form  of  four  distribution  functions 

curves  shown  in  Figure  1  have  the  same  mean  as  the  Weibull  curve.  The  Ex¬ 
tremal  Type  I  curve  of  Figure  1  was  derived  from  the  same  data  as  the  Weibull 
curve. 


Joint  Probability 


It  is  often  important  to  describe  an  event  by  more  than  one  variable, 
such  as  both  the  duration  and  peak  intensity,  in  which  case  the  joint  proba¬ 
bility  density  must  be  evaluated.  The  probability  that  variables  describing 
the  event  fall  within  specified  ranges  is  determined  from  the  joint  probabil¬ 
ity  density  in  a  similar  manner  as  with  single  variable  probability  density 
functions : 


P(*l  <  *,  <  I*;  .  «,  <  *2  <  K2  -•••.  V 


bl  h  b„ 


/  I  I  f(x1  ,  x2  . x^>  dXj 


dx....dx  (24) 
l  n 


a,  a 


1  “2  n 


and 


when  f(x.  ,  x»  . x  )  2:  0 

1  L  n 


OO  OO  00 


/  II  f(Xj  ,  x2  x^)  dx1  dx2...dxn  -  1 


(25) 


A  joint  distribution  function  can  be  defined  also: 


F(xt 


,  x  x  ) 

/  n 


■  f  f  f 


t  )  dt.  dt_...dt  (26) 
n  l  /  n 


The  marginal  probability  density  of  variable  x^  is  determined  by  integrating 
the  Joint  probability  density  function  over  the  entire  domain  of  all  variables 
except  x^  : 

00 

f(x1)  -  j  . . . 


00 

/ 


f  (x, 


x  )  dx, . . .dx.  . 
n  l  l—i 


dx 


x+1* 


.dx 


(27) 
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An  Important  feature  of  joint  probabilities  is  that  if  the  random  vari¬ 
ables  Involved  are  independent,  then  their  joint  distribution  function  is  the 
product  of  their  marginal  distribution  functions,  such  that: 

F(x.  ,  x  .  x  )  -  F(x.)  F(x  )...F(x  )  (28) 

i  L  n  l  i  n 

Another  Important  concept  of  joint  probabilities  is  conditional  proba¬ 
bility  density,  defined  in  the  case  of  two  random  variables  as  the  conditional 
probability  density  of  the  first,  x^  ,  given  that  the  second  takes  on  a  spec¬ 
ified  value,  x^  ,  or: 


f(x.  ,  x_) 

*1<X1  1  x2>  '  -  ’  lf  f(x2>  ■  0 


Conditional  distribution  functions,  such  as  F(x^  |  x^)  ,  also  can  be  defined, 
expressing  the  cumulative  probability  density  in  a  manner  analogous  to  single 
variable  density  functions.  Conditional  probability  densities  or  distribution 
functions  do  not  require  independence  for  their  definition. 

Concepts  Related  to  Evaluation  of  Risk 

A  traditional  measure  of  risk  of  encountering  an  event  of  a  specified 
intensity  x  ,  such  as  a  critical  flood  elevation,  wind  velocity,  or  wave 
height,  is  the  return  period,  RT(x)  .  This  is  defined  in  practical  terms  as 
the  average  waiting  period  between  exceedances  of  x  .  The  return  period  for 
variables  whose  rate  of  occurrence  is  Independent  of  their  intensity  (l.e., 
the  number  of  occurrences  per  unit  time  is  a  Poisson  process  with  a  mean  X) 
is  given  by  (Borgman  and  Reslo  1982): 


RT(x)  - 


U[1  -  F(x)  ]  } 


The  nonencounter  probability,  NE(x)  ,  is  defined  as  the  probability 
that,  during  a  specific,  time  Interval  L  ,  the  largest  intensity  encountered 
will  be  less  than  or  equal  to  x  .  This  can  be  expressed  In  terms  of  the 


distribution  function  F(x)  for  the  case  of  a  Poisson  process  as  (Borgman  and 
Reslo  1982): 


NE(x)  -  #~AL[ l-F(x) ]  (31) 

Expressed  In  terms  of  the  return  period: 

NE(x)  -  e“L/RT(x>  (32) 

This  last  relation  demonstrates  the  danger  of  misinterpreting  the  return 
period  as  a  frequency  of  occurrence  for  events  of  intensity  x  .  When 
L  -  RT(x)  ,  then  NE(x)  -  0.37  .  In  other  words,  there  is  a  63  percent  proba¬ 
bility  of  encountering  an  event  of  Intensity  x  during  the  time  Interval  L  . 
The  term  "risk"  Is  defined  as  the  probability  that  an  event  of  Intensity  x 
or  greater  will  occur  at  least  once  in  the  time  Interval  L  ,  which  Is 
1  -  NE(x)  . 

Another  concept  Important  In  risk  and  optimization  analyses  Is  that  of 
expectation,  E{x}  .  This  has  actually  already  been  defined  as  the  mean  of 
f(x)  : 


E{x)  -  p 


xtix)  dx 


(33) 


One  useful  feature  of  the  expectation  as  a  long-term  average  ol  the  values  of 
x  Is  that  the  expectation  of  a  function  of  x  ,  g(x)  can  be  defined  by: 


E[g(x) ) 


g (x)  f (x)  dx 


(34) 


Another  feature  with  respect  to  Poisson  processes  worth  noting  regards  the 
reference  time  period  for  risk  criteria,  such  as  estimation  of  the  average  an¬ 
nual  value  of  some  variable.  Relation  of  the  Poisson  parameter  X  to  expec¬ 
tations  of  functions  of  the  random  variable  x  (the  outcome  of  a  Poisson  pro¬ 
cess,  where  the  number  of  occurrences  per  unit  time  is  independent  of  the 
value  or  x)  is  easiest  demonstrated  by  an  example.  Assume  that  in  1  year  k 
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extreme  events  occur,  where  k  is  a  Poisson  variable.  Intensities  of  extreme 
events  are  represented  by  significant  wave  heights,  (1  -  1,  2,  3,  ...k)  . 

Damage  to  a  structure  caused  by  each  extreme  event  Is  assumed  to  be  a  function 
of  Hg  ,  D(Hg  )  .  Total  damage  In  the  year's  time  Is: 

k 

S?  ■  I  b(h.i>  (35) 

1-1 

Since  k  and  are  Independent,  then  the  expectation  with  respect  to 

Is: 

k 

E&)-  1  ’w.i”  (36> 

1-1 

Since  H  values  are  Independent  identically  distributed  random  variables, 

8  X 

they  all  have  the  same  expectation,  and: 

K??)  *  E(k)  EfD(H8)]  (37) 

Taking  the  expectation  of  k  to  be  the  average  number  of  extreme  events  per 
year  (-  the  Poisson  parameter,  X),  the  long-term  average  annual  storm  damage 
Is: 


E(f?)  '  1EtD(V! 


D(H  )f  (H  )  dH 

s  s  s 


(38) 


This  relation  Is  critical  In  optimization  ol  first  costs  against  estimates  of 
long-term  maintenance  costs. 

Regression  by  the  Method  of  Least  Squares 


An  Important  part  of  many  research  efforts  is  the  estimation  of  distri¬ 
bution  parameters  from  measured  data  by  regression  using  the  method  of  least 


squares.  Assumed  linear  relationships  between  an  Independent  variable  x  and 
a  dependent  variable  y  of  the  form: 

y  -  a  +  Bx  (39) 

can  be  tested  against  a  set  of  x  ,  y  data  and  the  differences,  e  ,  between 
the  estimated  y  and  the  predicted  value  measured.  These  differences  can  be 
due  to  measurement  errors  or  Inadequacies  In  the  assumed  relationship,  such  as 
neglect  of  other  Independent  variables  which  also  affect  the  value  of  y  . 

The  method  of  least  squares  allows  the  parameters  a  and  B  to  be  estimated 
by  constants  a  and  b  such  that  resulting  differences  In  the  predicted  ver¬ 
sus  measured  y  values  are  a  mlnumum.  Since  these  differences,  called  resid¬ 
uals,  could  be  both  positive  and  negative  and  therefore  have  a  tendency  to 
offset  each  other,  the  square  of  the  differences  is  minimized  instead.  Many 
nonlinear  relationships  can  be  transformed  Into  a  linear  form  to  take  advan¬ 
tage  of  this  technique. 

The  accuracy  or  reliability  of  least  squares  estimates  of  the  true  lin¬ 
ear  parameters  a  and  B  can  be  expressed  in  a  number  of  ways.  All  possible 

true  y  values  are  assumed  to  be  independently  normally  distributed  with 

2 

means  a  +  Bx  and  the  common  variance  o  .  Measured  values  then  can  be 
written  as: 


y1  -  a  +  Bxt  +  (40) 

where  e  represents  independent  normally  distributed  random  variables  with 
*  2 

zero  means  and  a  common  variance  a  .  This  variance  for  "n"  y  values  can 
be  estimated  in  terms  of  the  residuals  as: 


2 

8 

e 


I 


(a  +  bxi)  ]2 


(41) 


1-1 

where  s  Is  the  standard  error  of  estimate.  The  standard  error  is  in  units 
e 

of  y  and  represents  the  limit  within  which  approximately  68  percent  of  the 

absolute  values  of  all  errors  will  fall.  Another  quantitative  measure  of 

2 

variance  is  the  sum  of  the  square  residuals,  or  (n  -  2)  s 

e 


The  proportion  of  the  variation  of  y  values  which  can  be  attributed  to 
the  assumed  relationship  with  x  can  be  estimated  as  the  ratio  of  the  sum  of 
squared  residuals,  y  -  y  ,  to  the  sum  of  squared  deviations  of  y  from  the 
measured  mean,  y  ,  subtracted  from  1,  the  square  root  of  which  is  known  as 
the  nonlinear  correlation  coefficient,  r  : 


The  above  relation  has  the  advantage  over  other  correlation  formulas  that  it 
is  not  restricted  to  linear  relationships,  although  it  is  more  tedious  to 
compute . 

Confidence  that  can  be  placed  on  predictions  made  with  an  equation 
developed  by  the  least  squares  method  can  be  estimated  by  various  methods 
(Miller  and  Freund  1985,  Isaacson  and  MacRensle  1981).  The  upper  limit  of 
confidence  in  estimates  applied  as  design  criteria  always  should  be  addressed 
by  engineers  as  an  Integral  part  of  the  design  process,  particularly  if  pre¬ 
dictions  are  extrapolated  beyond  the  range  of  measured  data.  Techniques  for 
estimating  statistical  confidence  are  not  discussed  here  in  detail  since  this 
project  does  not  directly  involve  extrapolation.  It  should  be  noted,  however, 
that  obtaining  a  large  sample  Is  very  important  in  improving  statistical  con¬ 
fidence.  LeMehaute  and  Wang  (1984  and  1985)  have  made  special  note  of  the 
sensitive  effect  on  confidence  of  wave  statistics  attributable  to  the  number 
of  years  of  record  and  frequency  of  recordings.  Neglect  of  statistical  confi¬ 
dence  inherent  in  formulation  of  structural  design  criteria  can  lead  to  inade¬ 
quate  safety  and  higher  than  anticipated  maintenance  costs  for  structures  in¬ 
volved.  The  20  years  of  hlndcast  wave  data  at  3-hr  intervals  available  from 
the  WIS  program  are  valuable  in  this  regard. 


CHAPTER  III:  CHARACTERIZATION  OF  WAVE  CONDITIONS 


Basic  Sinusoidal  Concepts 

An  understanding  of  the  basic  theory  and  terminology  of  water  wave  me¬ 
chanics  is  necessary  for  interpretation  of  hlndcast  wave  information  and  any 
analytical  application  of  this  information.  Water  surface  waves  are  most  eas 
ily  described  as  wave  forms  of  sinusoidal  shape.  Certain  key  terms  with  ref¬ 
erence  to  this  simplified  concept  of  water  waves,  as  illustrated  in  Figure  2, 
include : 

1.  Wave  height,  H  -  the  vertical  distance  between  a  consecutive 
trough  and  crest 

2.  Wave  length,  L  -  the  horizontal  distance  between  two  consecutive 
crests  (or  troughs) 

3.  Wave  period,  T  -  visualizing  the  wave  form  as  travelling  horizon¬ 
tally,  the  time  for  two  consecutive  crests  (or  troughs)  to  pass  a  fixed  point 
usually  in  seconds 

4.  Wave  frequency,  f  -  nominally,  the  rate  at  which  consecutive 
crests  (or  troughs)  pass  a  fixed  point  (■  1/T),  in  hertz  (cycles  per  second) 


5.  Radial  frequency,  u>  -  the  radial  equivalent  of  frequency 
(u>  »  2ir/T) ,  also  in  hertz 

6.  Wave  number,  k  -  the  radial  equivalent  of  wave  length  (»  2x/L) 

7.  Phase,  $  -  the  radial  equivalent  of  the  horizontal  displacement, 

x'  ,  of  a  wave  crest  from  the  origin  of  the  reference  axis  at  time,  t  ■  0 
(-  2ir/x') 

The  basic  equation  which  defines  the  wave  profile  in  these  terms  is: 


n(x,  t)  -  ^  cos  Ckx  -  0t  +  ♦) 


where  n(x,  t)  is  the  lnstanteous  position  of  the  water  surface.  Considera¬ 
tion  of  the  sum  of  potential  and  kinetic  energy  inherent  in  a  travelling  wave 
of  this  form  (per  unit  surface  area)  can  be  estimated  by: 


where  p  is  the  mass  density  of  the  seawater.  This  total  energy  is  notably  a 
function  only  of  the  wave  height  squared  (Dean  and  Dalrymple  1984) . 

A  consideration  of  surface,  bottom,  and  transverse  boundary  conditions, 
with  simplifications  which  eliminate  all  but  first-order  differential  terms, 
yields  the  mathematical  equation,  known  as  the  dispersion  relation,  which  pre¬ 
dicts  effects  of  depth  on  wave  length: 


gk  tanh  (kd) 


where  g  is  the  acceleration  due  to  gravity  and  tanh  is  the  hyperbolic  tan¬ 


gent.  A  feature  of  sinusoidal  waves  which  is  consistent  with  this  relation  is 


2  2  2 

that  deepwater  wave  length,  Lo  =  (g/2’r)  T  **  5.12T  ft  or  1.56T  m.  The 


speed  at  which  a  wave  crest  travels,  the  phase  velocity,  C  ,  in  deep  water 


■  L  /T  «  5.12T  ft/sec  or  1.56T  m/sec.  The  change  that  occurs  in  shallower 
o 


water  is  that  wavelength  shortens  and  phase  velocity,  C  -  L/T  ,  increases. 
The  wave  height  also  is  affected,  first  slightly  decreasing,  then  increasing 
as  the  water  grows  more  shallow.  The  overall  tendency  of  water  waves  to 


change  form  as  depths  decrease  is  known  as  shoaling.  The  change  in  wave 
height  due  to  shoaling  is  governed  by: 


(46) 


where  H  and  H  are  shoaled  and  deepwater  wave  heights  and  K  is  the 
o  s 

shoaling  coefficient.  The  variable  C  is  the  shoaled  group  velocity,  the 

8 

speed  at  which  groups  of  waves  travel  which  is  also  the  speed  at  which  wave 
energy  approaches  shore: 

c8  ■  f  [•  -  sdrTjkd)]  <47> 

where  C  is  the  shoaled  phase  velocity  {=  L/T)  and  sinh  is  the  hyperbolic 
sin  function  (Dean  and  Dalrymple  1984). 

The  wave  form  becomes  steeper  in  decreasing  depths,  ultimately  reaching 
an  unstable  state  when  breaking  occurs.  The  point  at  which  breaking  actually 
occurs  is  not  fully  understood  at  this  time,  but,  based  on  the  theory  of  soli¬ 
tary  waves,  generally  occurs  at  the  point  where  the  wave  height,  H  =  0.78d  . 
Some  field  data  tend  to  show  that  most  locally  wind-generated  waves  (i.e. 
"seas")  break  in  deeper  water,  with  breaking  heights  on  the  order  of  0.6d  to 
0.7d  .  Very  long  waves  not  locally  generated  (i.e.  "swell")  may  not  break 
until  they  are  in  very  shallow  water,  however,  since  they  may  form  surging 
breakers  analogous  to  hydraulic  phenomena  known  as  "bores"  or  "hydraulic 
jumps. " 

The  discussion  above  is  meant  to  point  out  that  there  are  practical 
limits  to  wave  heights  at  most  coastal  sites  due  to  breaking,  but  that  these 
limits  are  as  yet  difficult  to  reliably  define  in  practice.  Furthermore,  sim¬ 
plifications  inherent  in  first-order  sinusoidal  theory  are  not  sufficiently 
accurate  for  engineering  purposes  in  many  shallow-water  situations  and  predic¬ 
tions  made  with  a  higher  order  wave  theory  must  be  applied. 

Shoaling  occurs  only  as  a  function  of  depth,  but  refraction  also  affects 
the  wave  form  as  a  function  of  wave  direction  with  respect  to  depth  contours 
of  the  sea  bottom.  Refraction  of  water  waves  is  analogous  to  refraction  in 
classical  physics  of  a  ray  of  light  passing  through  a  pane  of  glass  at  an 


angle.  The  most  frequently  observed  effect  of  water  wave  refraction  is  for 
waves  approaching  the  coast  at  an  angle  to  bend  around  as  their  crests  tend  to 
become  parallel  to  the  shoreline  in  shallow  water.  Snell's  Law  is  usually  ap¬ 
plied  to  describe  the  change  in  angle  of  water  waves  by  refraction  in  much  the 
same  way  as  it  is  in  optics,  commonly  stated  as: 


sin  0 


sin  0. 


(48) 


where  C  and  Cq  are  the  refracted  and  deepwater  phase  velocities  (■  L/T 
and  L^/T)  and  0  and  0q  are  the  refracted  and  deepwater  angles  of  wave 
crests  with  the  bottom  contours.  Snell's  Law  assumes  straight  and  parallel 
contours  between  deep  water  and  the  depth  at  which  the  above  relation  is  ap¬ 
plied.  The  relation  can  be  applied  in  increments  of  incident  versus  refracted 
angles  and  thus  applied  to  gently  curving  contours.  Refraction  usually  (ex¬ 
cept  in  cases  of  convergence  at  convex  contours)  causes  a  reduction  in  wave 
height,  which  is  superimposed  on  the  effect  of  shoaling,  according  to  the 
ratio: 


K  «  fCOB  V) 

r  H  \ cos  0  / 


1/2 


(49) 


where  H  and  Hq  are  the  refracted  and  deepwater  wave  heights  and  is 

the  refraction  coefficient  (Dean  and  Dalrymple  1984). 

Wave  diffraction  describes  the  effect  which  a  partial  barrier  has  on 
wave  heights  beyond  the  barrier.  It  is  the  process  which  allows  wave  energy 
to  leak  sideways  behind  an  obstruction  or  laterally  from  an  area  of  high  en¬ 
ergy  to  an  adjacent  area  of  lower  energy.  The  head  of  a  breakwater,  for  exam¬ 
ple,  will  cause  waves  to  diffract  behind  the  breakwater  into  its  geometric 
shadow,  even  though  it  may  prevent  any  other  form  of  wave  transmission. 

Larger  scale  landforms  and  submerged  formations  can  cause  a  degree  of  wave 
diffraction.  Precise  predictions  of  the  effects  of  diffraction  are  more  com¬ 
plicated  than  for  shoaling  and  refraction,  but  the  combined  effects  of  these 
three  forms  of  wave  transformation  are  important  in  explaining  observed  be¬ 
havior  of  water  waves  in  many  practical  situations.  The  complexity  of 


diffraction  often  requires  the  use  of  physical  scale  models  to  ensure  with 
confidence  satisfactory  performance  of  protective  structures  such  as  break¬ 
waters  enclosing  a  port  or  harbor  area. 


Irregular  Waves 

The  fact  that  real  ocean  waves  typically  appear  chaotic  with  little  reg¬ 
ular  form  was  mentioned  previously.  An  explanation  of  this  reality  is  that 
wave  groups  from  many  different  sources  with  different  heights,  periods, 
phases,  and  directions  are  interacting  in  the  small  area  we  observe  with  the 
resulting  superpositions  appearing  as  chaos.  Figure  3  illustrates  a 


hypothetical  point  in  time  when  two  sinusoidal  wave  groups  interact,  one  with 
50  percent  greater  height  and  period  and  a  ir/4  phase  difference.  The  waves 
would  appear  criss-crossed  when  viewed  from  above  if  their  directions  were  not 
parallel. 

Actually,  winds  that  create  the  waves  generate  a  range  of  heights  and 
periods.  Since  phase  velocity  varies  with  period,  longer  period  waves  travel 
faster  and  soon  leave  shorter  period  waves  behind.  Swell,  as  previously  de¬ 
fined,  refers  to  waves  which  have  completely  left  the  area  in  which  they  were 
generated.  These  waves  typically  have  periods  greater  than  about  9  or  10  sec, 
but  a  clear  distinction  does  not  exist.  Waves  which  are  still  within  the  in¬ 
fluence  of  the  generating  wind  system  are  called  "seas"  and  typically  are  dom¬ 
inated  by  shorter  period  waves  (less  than  9  sec) . 

The  distribution  of  individual  wave  heights  in  a  stationary  sea  state 
has  been  found  in  most  cases  to  follow  a  Rayleigh  distribution,  as  discussed 
in  the  previous  paragraphs  on  statistical  concepts.  Stationarity  technically 
is  the  condition  during  which  all  moments  (including  the  mean  and  variance) 


are  time  invariant  (Bendat  and  Piersol  1971).  A  small  sample  thus  can  be  ana¬ 
lyzed  and  taken  to  represent  the  entire  period  during  which  conditions  remain 
stationary.  Waves  at  sea  are  assumed  by  most  investigators  to  be  weakly  sta¬ 
tionary  for  periods  of  about  3  hr,  occasionally  for  as  much  as  6  hr,  but  sel¬ 
dom  longer.  This  is  more  of  a  tradition  related  to  the  practicalities  of  col¬ 
lecting  wave  data  than  a  precisely  defined  interval.  The  parameters  derived 
from  an  instanteous  measurement  (such  as  the  case  of  synoptic  hindcasting)  or 
from  a  20-min  recording  of  the  water  surface  elevations  are  therefore  typi¬ 
cally  taken  to  represent  a  much  longer  period  during  which  conditions  do  not 
change.  This,  of  course,  is  not  really  true,  but  as  long  as  the  changes  are 
not  drastic  and  are  generally  within  the  confidence  limits  of  the  statistical 
parameters  of  interest,  this  practice  is  acceptable. 

Wave  periods  do  not  lend  themselves  as  readily  as  do  wave  heights  to 

representation  by  a  standard  statistical  distribution  such  as  the  Raleigh 

distribution.  Bretschneider  (1959),  however,  found  that  the  distribution  of 
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squared  wave  periods,  T  ,  for  seas  followed  a  Rayleigh  distribution.  Other 
investigators  have  applied  a  variety  of  standard  distributions,  and  special¬ 
ized  empirical  distributions  also  have  been  developed. 

The  practice  of  coastal  engineers  in  the  last  10  years  has  largely 
shifted  from  considerations  of  wave  period  exclusively  in  the  time  domain  to 
frequency  domain  considerations.  Decomposition  of  a  time  series  of  water  sur¬ 
face  elevations  into  a  set  of  Incremental  sinusoids,  each  represented  by  an 
amplitude  (=  H/2)  and  a  frequency  (-  l/T),  can  be  accomplished  by  transforma¬ 
tion  of  the  time  series  into  its  equivalent  Fourier  series.  Wave  conditions 
thus  can  be  represented  by  the  distribution  of  wave  energy  (proportional  to 
amplitude  squared  per  Equation  44)  as  a  function  of  frequency,  or  a  wave 
spectrum. 

Figure  4  illustrates  a  wave  spectrum  with  two  "peaks,"  one  representing 
swell-type  waves  and  the  other  representing  coexistent  seas.  The  inverse  fre¬ 
quency  of  the  dominant  peak  is  in  practice  usually  taken  as  the  peak  period, 
which  is  generally  assumed  as  the  most  probable  period  in  the  sea  state.  This 
is  a  "one-dimensional"  spectrum  which  does  not  account  for  the  direction  of 
wave  energy  propagation.  More  complex  procedures  have  been  developed  to  ex¬ 
press  the  distribution  of  wave  energy  as  a  function  of  both  frequency  and  di¬ 
rection.  The  most  common  practice  is  to  treat  the  directional  spread  of  wave 
energy  to  be  independent  of  the  distribution  of  energy  by  frequency.  This 
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Figure  4.  An  example  of  a  double-peaked 
energy  density  spectrum 


allows  application  of  a  spreading  function  0(0)  which,  when  multiplied  by 
the  one-dimensional  spectrum  S(f)  ,  yields  the  two-dimensional  spectrum 
Sff,  6)  : 


S(f,  6)  -  S(f)  0(0) 


(50) 


The  form  of  a  spectrum  Is  quite  sensitive  to  the  analytical  procedures 
applied,  particularly  "smoothing"  performed  to  Improve  statistical  confidence 
at  the  cost  of  resolution.  Most  spectral  analysis  procedures  actually  deal 
with  discrete  frequencies  (**  2u/T  of  the  individual  sinusoids)  which,  when 
averaged  over  equal  intervals,  yield  a  smoother  looking  plot  with  more  narrow 
confidence  bands.  A  jagged  looking  spectrum  will  have  wider  confidence  limits 
than  a  smoothed  spectrum  computed  from  the  same  data. 


Integration  of  a  wave  spectrum  which  has  been  computed  as  energy  per 
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frequency  band,  E/Af  (e.g.  m  / Hz)  ,  versus  frequency  yields  the  total  energy 
of  the  sea  state.  This  relates  directly  to  actual  variance  of  the  water  sur¬ 
face  elevations  such  that: 
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vs 


/■ 


S(f)  df 


(51) 


where  o  is  the  variance  of  the  water  surface  elevations  and  S(f)  is  the 
ws 

computed  energy  density  spectrum.  Spectra  in  this  form  are  often  taken  as 
continuous  functions  since  it  is  reasonable  to  expect  wave  energy  to  be  gener¬ 
ated  in  continuous  frequencies. 

A  parameter  in  units  of  wave  height  which  has  been  used  to  represent  the 
range  of  wave  heights  in  a  sea  state  is  the  zero  moment  wave  height, 

H  ■  4o  .  The  "zero  moment”  title  comes  from  integration  of  fnS(f)  with 
respect  to  f  where  n  ,  the  power  of  f  in  the  Integral,  is  zero  as  with 
Equation  51.  This  wave  height  has  been  found  to  be  very  close  to  the  signifi¬ 
cant  wave  height,  Hg  ,  of  Rayleigh  distributed  seas  in  deep  water.  Hg  typi¬ 
cally  departs  from  H  in  shallow  water  (Thompson  and  Vincent  1983).  The 

mo 

zero  moment  wave  heights  corresponding  to  two  interacting  wave  groups  of 
double-peaked  energy  density  spectra,  as  illustrated  in  Figure  4,  can  be 
estimated  by  splitting  the  spectrum  between  peaks  and  integrating  each  side 
separately.  There  is  no  widely  accepted  way  to  estimate  the  parameters  of 
multiple  wave  groups  from  their  combined  spectrum,  but  this  method  gives  an 
indication  of  their  relative  intensity  as  potential  structural  design 
criteria. 

A  number  of  parameterized  spectra  have  been  developed  in  the  effort  to 
rebate  wave  conditions  to  winds  and  geographical  factors  which  constrain  gen¬ 
eration  of  waves  at  sea.  These  parametric  spectral  forms  nearly  all  apply  to 
waves  in  the  generation  phase,  i.e.  seas,  not  swell.  The  four  most  important 
factors  in  wave  generation  are  wind  velocity  (and  resultant  stress)  over 
water,  duration  of  that  velocity,  fetch  (distance  over  water  which  the  wind 
blows),  and  water  depth.  Depth  limitations  on  wave  spectra  are  the  most 
recent  effects  to  be  reliably  defined  in  combination  with  other  primary  con¬ 
straints.  Other  factors  which  also  can  be  significant  are  preexisting  waves 
(wave-wave  interaction)  and  the  presence  of  strong  currents  (wave-current 
interaction).  Waves  generated  by  winds  of  a  given  velocity  in  water  of  a 
given  depth  thus  are  either  duration  limited,  fetch  limited,  or  fully  devel¬ 
oped  and  may  be  affected  by  waves  coming  into  the  generation  area  from  a  dis¬ 
tant  source  and  strong  currents.  Virtually  all  parametric  spectral  shapes 
have  the  "tail"  of  the  spectrum,  the  portion  to  the  right  of  the  peak, 
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proportional  to  f  ^  ,  following  the  work  of  Phillips  (1977).  An  advanced 
form,  as  an  example,  is  the  TMA  spectrum  (Hughes  1984),  which  includes  the 
depth  limitation: 

7  «■  ,  — 5/4(f/f  )"4  exp  -(f/f  -l)2/2 al 

S(f,  d)  -  og  r-3  (2ir)"%(2irf,  d)e  P  Y  P  (52) 

where  <fr(2irf,  d)  is  a  function  of  depth  (d)  ,  k  (the  wave  number,  2it/L)  , 
and  u)  (the  radial  frequency,  2ir/T)  allowing  portions  of  the  spectrum  to  be 
transformed  by  linear  wave  theory.  The  term  a  is  the  Phillips  equilibrium 
constant,  which  has  recently  been  taken  to  be  a  function  of  depth,  wind  speed, 
and  peak  frequency,  f  .  The  Y  term  is  the  "shape  parameter"  which  is  a 
function  of  wind  speed  and  fetch.  The  term  is  an  empirical  factor  af¬ 

fecting  shape  of  the  spectrum  on  either  side  of  the  peak.  This  form  applies 
to  fully  developed  or  "saturated"  seas  in  decreasing  depths.  Figure  5  illus¬ 
trates  the  effect  of  changing  depth  on  TMA  spectral  shape.  The  deepwater  pre¬ 
decessor  of  the  TMA  spectrum,  the  JONSWAP  spectrum,  now  is  widely  used  to  pre¬ 
dict  both  fetch  and  duration  limited  wave  growth  in  deep  water  (Vincent  1984). 


CHAPTER  IV:  WAVE  INFORMATION  STUDIES  HINDCAST  DATABASF 


General  Background  of  Phases  I  and  II 

The  WIS  program  of  the  US  Army  Engineer  Waterways  Experiment  Station 
began  In  1976  with  the  goal  of  providing  a  long-term  (20-vear)  hlndcast  of 
wave  Information  for  use  In  development  of  design  criteria  for  coastal  proj¬ 
ects.  The  term  "hlndcast"  refers  to  the  technique  of  simulating  historical 
wind  and  wave  generation  from  pressure  data  available  from  surface  weather 
charts.  The  basic  raw  data  for  hindcasting  thus  are  lnstanteous  pressure  re¬ 
cordings  which  meteorologists  have  applied  to  produce  pressure  fields  delin¬ 
eated  by  Isobars  and  other  notation  common  to  surface  weather  charts.  These 
"highs,"  "lows,"  "fronts,"  "troughs,"  and  "ridges"  are  then  applied  to  simu¬ 
late  the  effect  of  corresponding  wind  fields  on  the  surface  of  the  ocean. 

The  WTS  program  first  transcribed  Into  digital  form  pressures  from  sur¬ 
face  weather  charts  from  1956-1975  for  the  North  Atlantic,  Gult  of  Mexico,  and 
North  Pacific,  with  as  much  checking  for  accuracy  and  consistency  as  the  basic 
data  allowed  (Corson,  Resio,  and  Vincent  1980).  This  information  was  avail¬ 
able  at  6-hr  intervals.  Winds  which  would  have  existed  with  each  consecutive 
distribution  of  pressures  next  were  simulated  by  a  series  of  numerical  models 
assuming  quasigeostrophic  flows  and  a  planetary  boundary  layer  which  yielded 
surface  level  (19.5-m  elevation)  wind  fields.  These  wind  fields  were  in  turn 
adjusted  with  observations  of  actual  wind  velocities,  wherever  possible 
(Resio,  Vincent,  and  Corson  1982). 

Given  the  database  of  surface  level  winds  created  by  the  steps  above, 
basin  geometry  and  grid  were  defined  for  numerical  simulation  of  deepwater 
wave  generation.  Figures  A-l  and  A-2  illustrate  deepwater  (Phase  I)  grids  for 
the  North  Atlantic  and  North  Pacific  Oceans.  Execution  of  a  deepwater  numeri¬ 
cal  model  of  wave  generation,  which  took  into  account  fetch,  duration,  direc¬ 
tional  spreading  effects,  and  wave-wave  interaction,  produced  a  database  of 
two-dimensional  spectra  and  related  parameters  at  intersections  of  the  grid 
lines.  Detailed  wave  information  was  retained  only  at  intersections  marked 
with  dots  and  published  in  written  form  (Corson  et  al.  1981  and  Ragsdale  1983) 
for  the  numbered  sites. 

Phase  II  of  the  WTS  program  performed  simulations  at  3-hr  intervals  of 
wave  generation  in  a  manner  similar  to  Phase  I  (deep  water),  but  at  a  finer 


scale  and  In  transitional  depths  of  the  continental  shelf.  Figures  A-3  and 
A-4  show  the  Atlantic  and  Pacific  Phase  ll  grids  and  stations  where  wave  In¬ 
formation  has  been  published  (Corson  et  al.  1982  and  Ragsdale  1983).  In  ad¬ 
dition  to  Phase  I  factors.  Phase  II  simulations  took  Into  account  the  shelter¬ 
ing  effect  of  large-scale  land  masses,  refraction,  and  shoaling.  The  Phase  I 
wave  Information  served  as  a  boundary  condition  at  the  seaward  limit  of  the 
Phase  II  grid. 

Neither  Phase  I  nor  Phase  II  distinguished  seas  and  swell,  but  rather 
dealt  with  Individual  discrete  frequency  bands  over  the  entire  two-dimensional 
spectrum  at  any  point.  Phase  III  decomposed  this  spectrum  into  seas  and 
swell,  treating  seas  as  two-dimensional  spectra  and  swell  as  monochromatic, 
unidirectional  wave  groups.  The  definition  of  swell  as  waves  which  have  trav¬ 
elled  beyond  the  area  in  which  they  were  generated  was  applied.  This  approach 
economized  computations  by  taking  advantage  of  the  fact  that  swell  typically 
has  its  energy  highly  concentrated  In  a  narrow  band  of  frequencies,  which  is 
close  to  a  monochromatic  condition.  Wave  parameters  computed  and  recorded  In 
the  Phase  III  database  Included  zero  moment  wave  height,  peak  period,  and  dom¬ 
inant  direction  of  propagation.  Monochromatic  equivalents  were  recorded  In 
the  case  of  swell  and  combined  wave  heights  were  recorded  as: 

combined 

Period  and  direction  recorded  In  the  "combined"  category  corresponded  to  the 
peak  period  and  dominant  direction  of  either  the  sea  or  swell,  whichever  had 
the  higher  zero  moment  wave  height  (Brooks  and  Corson  1984).  The  Phase  III 
approach  is  most  valid  for  coasts  with  straight  and  parallel  contours  and  is 
leas  precise  in  more  complex  bathymetry. 

Phase  III  Shallow-Water  Wave  Information 

Phase  III  efforts  of  the  WIS  program  were  directed  at  providing  wave  in¬ 
formation  suitable  as  design  criteria  for  a  great  many  coastal  endeavors  In  a 
depth  of  10  m  at  10-mile  (16.1-km)  Intervals  along  the  Atlantic  (Jensen  1983a) 
and  Pacific  (Ragsdale  1983)  coasts  of  the  continental  United  States.  This 
task  dealt  with  transformation  of  wave  conditions  from  Phase  II  stations  to 


166  Atlantic  and  134  Pacific  Phase  III  stations.  Figure  A-5  illustrates  a 
section  of  the  Atlantic  Phase  III  stationing  system  and  adjacent  Phase  II  sta¬ 
tions.  The  magnitude  of  data  processing  requirements  and  complexity  of  the 
coast  at  this  finer  scale  led  to  procedures  for  estimating  wave  conditions  in 
shallow  water  (10  m  depth)  described  briefly  below. 

A  spectral  (frequency  domain)  approach  to  wave  transformation  was  sought 
to  reduce  computational  time  required  to  simulate  wave  transformation  in  the 
time  domain.  A  parameterized  spectrum  was  necessary  for  this,  but  one  as  com¬ 
plex  as  the  TMA  spectrum,  or  the  most  refined  spectral  forms  available  at  the 
time  of  the  Phase  III  procedure  formulation,  would  have  provided  an  unmanage¬ 
able  computational  burden.  The  one-dimensional  parameterized  spectrum  chosen 
for  Phase  Ill  simulations  had  the  following  form: 


2  -5  -4 

S(f )  -  «g  f  3(2*) 


for  f  *  f 


(54) 


S (f )  ■  ag2fp5(2ir)  4  exp 


f'-ftf 


for  f  <  f 


(55) 


which  applied  the  well-accepted  f  right-hand  tail,  but  limited  free  param¬ 
eter  determination  to  only  two  variables,  a  and  f  (Kitaigordskii  1962). 

P 

A  spreading  function,  assumed  to  be  independent  of  the  one-dimensional 
spectral  form,  was  defined  as: 


0(0)  -  I-  cos4  (9-6') 

J7T 


(56) 


where  0'  is  the  predominant  direction  of  propagation.  Thus,  the  two- 
dimensional  form  was: 


S(f,  0)  -  S(f)O(0) 


(57) 


Within  each  10-mile  (16.1-km)  interval  defined  as  Phase  III  stations 
along  the  coast,  bottom  contours  were  assumed  to  be  straight  and  parallel.  A 
specific  orientation  was  assigned  to  each  interval  such  that  departure  of  this 
assumption  from  the  true  situation  was  minimized.  The  processes  of  refraction 


and  shoaling,  as  defined  by  Snell's  Law  and  sinusoidal  theory,  were  applied  to 
increments  ot  frequency  and  direction  of  the  directional  distribution  defined 
by  S(f,  0)  .  Wave  energy  propagating  seaward  was  ignored. 

The  geometric  relationship  between  a  Phase  III  station  and  adjacent 
Phase  II  stations  from  which  the  model  derived  its  input  was  the  most  impor¬ 
tant  consideration  in  addressing  sheltering  in  Phase  III,  Basically,  the  geo¬ 
metric  shadow  of  a  landform  to  wave  energy  from  a  specific  direction  was  con¬ 
sidered  as  absolute,  i.e.,  no  energy  was  propagated  into  the  shadow  area. 

This  is  a  gross  simplification,  but  it  made  the  simulation  of  sheltering  ef¬ 
fects  practical  for  Phase  ill.  Discrete  combinations  of  frequency  and  direc¬ 
tion  were  considered  incrementally  with  respect  to  sheltering,  as  they  were 
with  refraction  and  shoaling. 

The  problem  of  wave-wave  interaction  and  the  losses  it  can  cause,  evi¬ 
denced  by  white  caps  and  other  signs  of  turbulent  energy  dissipation,  was  ad¬ 
dressed  by  definition  of  another  spectral  form  for  shallow  water.  Principles 
of  similarity  were  applied  to  derive  a  form  consistent  with  Phase  I  and  II 
deepwater  considerations,  which  predicted  the  spectrum  in  shallow  Water: 


S(f)  -  agh(8ir)"2f"3 


for  f  £  f 


This  relation  is  consistent  with  the  visualization  that  energy  losses  due  to 
wave-wave  interaction  tend  to  occur  at  high  frequencies,  while  energy  at  lower 
frequencies  is  conserved.  A  further  application  of  equilibrium  principles 
allowed  derivation  of  an  integrated  form  of  this  equation  which  describes  the 
dependency  of  sea  wave  heights  on  depth: 
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where  f  ■  0.9f  is  a  energy  cutoff  frequency  (lower  integration  limit)  and 
c  p 

(H  )  is  the  upper  limit  of  seas  wave  heights.  Surf  zone  breaking  was 
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max 

treated  differently  for  swell,  however,  in  the  manner  of  estimating  breaker 
heights  for  monochromatic  waves.  A  breaking  coefficient  of  0.6  was  applied, 
which  is  consistent  with  recent  measurements  of  breaking  waves  by  the  WES 
(Jensen,  Robert  E. ,  verbal  communication,  February  1986): 


0.6d 


(60) 


H 

max 

Extensive  comparisons  have  been  made  between  the  limited  measured  wave 
data  available  and  WIS  wave  information,  generally  with  acceptable  results 
(Corson  and  Resio  1981).  The  reduction  of  measurements  made  by  wave  gages 
also  involves  compounded  assumptions ,  and  discrepancies  between  wave  informa¬ 
tion  based  on  gage  data  and  Phase  III  wave  information  could  not  always  be  re¬ 
solved.  More  accurate  techniques  are  available  for  site-specitic  simulation 
of  the  transformation  of  waves  into  shallow  water.  These  methods  unfortu¬ 
nately  were  too  complex  to  apply  systematically  on  the  scale  of  the  WIS 
Phase  III  endeavor,  though  improvements  are  under  consideration.  The  pres¬ 
ently  available  end  product  of  Phase  TII  is,  however,  an  excellent  tool  for 
coastal  engineers  to  use  in  the  planning  and  preliminary  design  stages  of 
coastal  projects  for  development  of  design  criteria.  More  complex  and  expen¬ 
sive  numerical  simulations  and  physical  scale  models  can  be  performed  in  the 
detailed  design  phase  after  the  economic  feasibility  and  f inanceability  of  the 
project  has  been  ensured.  Even  in  the  final  stage,  some  basis  of  experiment 
design  and  cross-check  on  other  sources  of  wave  information  is  necessary.  The 
20-year  period  of  record  for  the  WIS  database  can  rarely  be  exceeded  by  other 
reliable  sources.  The  WIS  wave  information  provides,  therefore,  a  vast  im¬ 
provement  to  the  confidence  of  each  design  effort  to  which  it  is  applied. 
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CHAPTER  V:  LITERATURE  REVIEW  OF  STORM  DURATION  STUDIES 


Recent  Literature  on  the  Duration  of  Sea  States 

Table  1  presents  mean  durations  for  various  weather  types  in  the  British 
isles  which  were  excerpted  from  Barry  and  Perry  (1973).  The  weather  type 


Table  1 

Mean  Durations  of  Weather  Types  in  the  British  Isles 


Mean  Duration 

(days) 

January 

July 

Weather  Type 

1910-1930 

1948-1968 

1910-1930 

1948-1968 

Westerly 

4.1 

1.7 

2.6 

2.7 

Northerly 

1.5 

1.8 

2.0 

1.8 

Easterly 

1.9 

2.6 

2.0 

1.8 

Southerly 

2.0 

1.8 

1.7 

1.3 

Cyclonic 

1.4 

1.3 

1.7 

1.9 

Anticyclonic 

2.2 

1.9 

2.5 

2.2 

identified  as  "cyclonic"  is  assumed  to  meet  the  standard  definition  of  winds 
circulating  around  a  low  pressure  area  (Lester  1973) ,  corresponding  to  the 
extratropical  cyclonic  events  which  are  simulated  in  the  WIS  program.  This 
type  of  weather  is  noted  to  have  a  mean  duration  of  1  to  2  days  in  Great 
Britain,  with  some  seasonal  variation.  Statistics  of  this  type  would  surely 
vary  from  region  to  region,  but  the  order  of  magnitude  in  hours,  say  less  than 
100  but  more  than  10,  can  serve  in  this  investigation  of  storm  characteristics 
as  a  rough  first  measure  of  a  reasonable  mean  duration.  The  untrained  intui¬ 
tion  of  any  regular  viewer  of  television  weather  reports  would  likely  agree 
with  this  typical  range. 

Surprisingly  little  material  was  available  in  the  coastal  engineering 
and  oceanographic  literature  which  dealt  directly  with  the  duration  of  extreme 
events  at  sea  or  of  extreme  wave  conditions.  Occasional  references  were  made 
to  a  3-hr  period  of  wave  height  stationarity  assumed  for  practical  purposes  in 
measurement  programs  (e.g.,  Agerschou  et  al.  1983  and  Massie  1976).  The  in¬ 
terval  between  samples  of  wave  measurements  is  commonly  set  at  3  hr. 


Publications  of  WIS  wave  information  (Corson  et  al.  1981  and  1982  and  Jensen 
1983a)  tabulated  durations  of  significant  wave  heights  above  selected  thresh¬ 
olds,  but  did  not  discuss  trends  or  other  implications  inherent  in  this 
information. 

North  Sea  Investigations  of  Houmb  and  Vik 

The  most  rigorous  work  to  date  has  been  a  series  of  studies  by  two 
Norwegian  investigators  (Houmb  1971,  Houmb  and  Vik  1975,  Vik  and  Houmb  1976, 
and  Houmb  and  Vik  1977).  Other  authors  have  reviewed  this  work  (e.g.  Battjes 
1977,  PIANC  1979,  and  Bruun  1985),  but  no  significant  advances  seem  to  have 
been  made  regarding  the  characterization  ot  extreme  event  durations  following 
Houmb  and  Vik  (1977).  Their  work  on  the  duration  of  sea  states  culminated  in 
the  findings  of  the  last  reference,  which  will  ba  reviewed  in  detail  in  the 
following  paragraphs. 

Houmb  and  Vik  (1977)  considered  both  the  duration  of  extreme  events, 
specified  as  the  time  during  which  the  significant  wave  height  exceeded  a 
given  threshold,  and  the  duration  of  "calms"  between  these  extremes.  The 
basis  of  their  investigations  was  wave  recordings  made  at  five  North  Sea  sites 
where  depths  varied  from  80  to  250  m.  Three  sites  involved  time  series  mea¬ 
surements  made  for  20  min  every  3  hr.  A  fourth  site  involved  10-min  time 
series  measured  every  4  hr.  The  sequences  of  these  measurements  were  not  con¬ 
tinuous  and  varied  in  total  period  of  record  from  3  to  31  months.  The  fifth 
site  provided  observations  from  a  rescue  vessel  every  3  hr  from  1959  to  1974 
during  October  through  March  only.  These  observations  classified  predominant 
wave  heights  into  classes  of  0.5  m. 

A  theoretical  approach  toward  prediction  of  variation  of  storm  durations 
was  first  proposed  by  Houmb  and  Vik  (1977)  which  took  the  frequency,  or 
marginal  probability  density,  of  threshold  up-crossings  (i.e.  H'  *  dH  /dt 

8  S 

was  positive)  as: 

f (H  )  -  H ' f (H  ,  H')  dH'  (61) 

t  S  t  S  8 


where  H  is  the  specified  threshold  and  f(H  ,  H')  is  the  joint  probability 
t  t  8 

density  of  Hg  and  its  time  derivation,  H^  .  The  average  duration  of  ex¬ 
treme  events,  t(H  )  ,  (H  >  H  )  was  derived  to  be: 
t  St 


(62) 


L[1  -  F(Ht>]  [1  -  F(Ht)] 

t(Ht)  fOT)I  fClT) 

where  L  is  the  period  of  interest  (say  50  years)  and  F(Ht)  is  the  cumula¬ 
tive  distribution  of  H  evaluated  at  H  ,  or  the  probability  that  H  is 

s  t  s 

equal  to  or  less  than  .  The  quantity  [1  -  F(Ht)]  is  the  probability  that 

H  is  greater  than  H  .  The  average  number  of  up-crossings,  i.e.  the  aver- 
s  t 

age  number  of  extreme  events,  during  the  period  L  was  taken  to  be  f(Ht)L  , 
where  f(H  )  is  the  probability  density  of  H  at  H  given  above. 

t  S  t 

The  rate  at  which  Hg  changes  (from  one  stationary  period  to  the  next) 

was  assumed  to  be  a  Poisson  process,  i.e.  H'  was  assumed  to  be  independent 

s 

of  .  The  joint  probability  density  function  f(Hg,  H^)  could  then  be 

evaluated  as: 


f (H  ,  H’)  ■=  f(H  )f(H')  (63) 

s  s  s  s 

The  marginal  probability  density  function  f(Hg)  was  assumed  to  follow 
a  Weibull  distribution  whose  corresponding  distribution  function  had  the  form: 


and 


(64) 


T(H  H  )T-1 

f  (H  )  -  - — - f  exp 

8  (H  -  H  ) 
c  o 


(65) 


where  H  ,  H  ,  and  T  are  parameters  of  the  distribution, 
c  o 

The  function  f(H')  was  assumed  to  be  normally  distributed  with  zero 
s 

mean  for  positive  values  of  H'  (increasing  H  ).  The  data  seemed  to  support 

S  8 

this  assumption.  This  gave  f(H^)  as: 


f(H;) 


(66) 
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where  o.  is  the  standard  deviation  of  H'  which  was  evaluated  from  the 
h  s 

data.  The  differences  between  values  computed  for  increasing  and  de¬ 

creasing  H  were  found  to  be  negligible.  Furthermore  o,  was  not  noted  to 
s  n 

follow  a  seasonal  pattern.  This  application  of  the  above  normal  distribution 
with  zero  mean  gave  the  advantage  of  requiring  only  one  parameter,  ,  to  be 


determined  empirically,  in  addition  to  those  (H 


H  ,  and  T)  for  F(H  )  . 
o  8 


The  resulting  function  for  the  mean  duration  t(H^)  reduced  to: 


t  (H  ) 


2tt(H  -  H  ) 
c  o 

To.  Or  -  h  ) 
h  t  o 


The  cumulative  distribution  of  measured  durations  was  found  to  be  well 
represented  by  a  Weibull  distribution  of  the  form: 


F(t)  =  1  -  exp 


where  a  is  the  shape  parameter  and  t  is  the  scale  parameter.  Average 

c 

durations  estimated  by  the  t(Ht>  function  derived  above  also  compared  well 
with  means  computed  from  the  set  of  measured  durations.  Houmb  and  Vik  (1977) 
gave  examples  of  how  this  formulation  could  be  applied  in  the  conduct  of  off¬ 
shore  oil  explorations,  as  in  prediction  of  duration  of  operation  down  time 
caused  by  extreme  wave  conditions. 

The  formulation  of  Houmb  and  Vik  (1977)  was  well  defended  in  terms  of 
conceptual  limits  ot  parameters  such  as  and  .  They  tested  their 

hypotheses  as  well  as  possible  with  their  limited  data  set,  but  urged  in  their 
conclusions  that  further  investigations  be  pursued  with  more  comprehensive 
wave  information. 


CHAPTER  VI:  EXTREME  EVENT  IDENTIFICATION 


Choice  of  Sites 


Each  Phase  III  site  includes  58,440  records  of  wave  information  3  hr 
apart  from  0000  (midnight)  January  1,  1956,  to  2400  (midnight)  December  31, 

I  1975  (20  years).  Four  sites  were  originally  chosen  for  analysis,  two  on  the 

'  Atlantic  coast  and  two  on  the  Pacific  coast.  A  third  Atlantic  site  was  later 

chosen  when  it  was  discovered  the  first  two  had  very  similar  distributions  of 
significant  wave  heights.  The  five  sites  ultimately  investigated  are  listed 
in  Table  2.  They  were  Intended  to  represent  a  wide  geographical  spread  in 

Table  2 

i 

!  WIS  Phase  III  Stations  Investigated 


Station 

Site 

Latitude 

Longitude 

A3061 

Atlantic  City,  New  Jersey 

39.34°  N 

74.47°  W 

A3083 

Nagshead,  North  Carolina 

35.94°  N 

75.61°  W 

A3142 

Daytona  Beach,  Florida 

29.20°  N 

81.00“  W 

P3036 

Newport,  Oregon 

43.63°  N 

124.08“  W 

P3105 

Half-Moon  Bay,  California 

37.45°  N 

122.45°  W 

hopes  that  analysis  would  reveal  any  important  universal  traits  or  significant 
geographical  differences.  Figure  6  shows  their  relative  location  along  the  US 
coasts.  Statistics  published  by  the  WIS  program  (Jensen  1983a,  b)  for  the 
Atlantic  sites  are  presented  in  Appendix  A.  Wave  height  frequency  tables  (not 
yet  published  by  the  WIS  program)  for  the  two  Pacific  sites  also  are  presented 
in  Appendix  A. 


Basic  Treatment  of  WIS  Phase  III  Wave  Information 


Table  A1  illustrates  format  and  unit  conventions  of  the  WIS  Phase  III 
database.  Dates  are  given  as  year/month/day  and  times  referenced  to  the  24-hr 
clock  (i.e.,  military  time).  Wave  heights,  i.e.  the  zero  moment  wave  heights 
derived  for  each  3-hr  time  step,  are  reported  in  centimetres.  Wave  periods. 
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Figure  6.  Geographic  relation  of  sites  investigated 

i.e.  the  peak  periods  of  the  hindcast  spectra,  are  reported  to  the  nearest 
second.  Direction  or  azimuth  is  reported  in  degrees  relative  to  the  shore¬ 
line,  such  that  90  deg  is  a  wave  direction  travelling  straight  into  the 
straight  and  parallel  contours  assumed  for  each  10-mile  (10.1-km)  shoreline 
increment.  Combined  statistics  presented  in  Table  A1  include  the  geometric 
average  wave  height  (Equation  53)  and  the  peak  period  and  predominant  direc¬ 
tion  of  either  seas  or  swell,  whichever  had  the  highest  zero  moment  wave 
height.  Combined  statistics  were  applied  in  analyses  of  this  study,  though 
they  were  not  actually  a  part  of  stored  wave  information  and  had  to  be  com¬ 
puted.  Mean  and  maximum  duration  of  exceedance  of  selected  wave  heights  were 
reported  for  the  Atlantic  sites  by  Jensen  (1983b)  and  are  included  in  Appen¬ 
dix  B.  A  comparison  of  those  statistics  with  results  from  this  investigation 
is  made  later  in  this  report. 

The  Problem  of  Extreme  Event  Identification 

The  work  of  Houmb  and  Vik  (1977)  on  duration  of  sea  states  was  appar¬ 
ently  performed  exclusively  with  significant  wave  heights  crossing  an 


arbitrary  threshold.  This  implies  that  the  significant  (or  zero  moment)  wave 
height  is  the  most  appropriate  measure  of  the  extreme  events'  intensities  for 
applications  of  duration  statistics.  Other  parameters  can  be  conceived,  how¬ 
ever,  which  might  be  better  representatives  of  the  overall  intensity  or 
extreme  nature  of  a  storm.  The  most  obvious  alternate  parameter  would  be  peak 
period,  to  which  refraction,  shoaling,  and  wave  breaking  are  all  quite  sensi¬ 
tive.  Wave  length  might  be  another,  although  wave  length  at  any  depth  is  a 
function  of  period.  Wave  steepness,  H/L  ,  is  commonly  associated  with 
breaker  characteristics  and  forces  on  coastal  structures.  If  the  ratio  of 

zero  moment  wave  height  to  deepwater  wavelength  corresponding  to  the  peak  pe- 

2 

riod  is  used,  representative  wave  steepness  becomes  2irH/gT  .  The  2tt  fac¬ 
tor  is  commonly  dropped  as  a  part  of  this  dimensionless  steepness  parameter  in 
favor  of  H/gT^  . 

2 

Wave  severity,  H  L  ,  has  recently  become  of  interest  as  a  factor 
closely  related  to  stability  of  rubble-mound  structures  (Graveson  et  al.  1980 
and  Ahrens  1984).  Wave  severity  can  be  thought  of  as  the  ratio  of  wave  height 
cubed  (the  traditional  wave  parameter  for  evaluation  of  rubble-mound  stabil¬ 
ity)  to  wave  steepness,  H/L  .  Again,  significant  or  zero  moment  wave  height 

and  deepwater  wave  length  corresponding  to  the  peak  period  of  the  spectrum  are 

2  2  2 

used  for  convenience,  yielding  H  L  ■  2irH  / gT  .  It  should  be  noted  that  the 

four  parameters  discussed  so  far  vary  the  relative  influence  of  wave  height 

2  2  2 

and  period  in  the  following  order:  H  ,  T  ,  H/T  ,  and  H  /T  .  These 
parameters  also  could  be  used  to  define  extreme  event  duration  as  the  time 
during  which  consecutive  parameter  values  exceed  a  specified  threshold  value. 

A  fifth  parameter  which  might  be  important  with  respect  to  duration  of 
extreme  wave  conditions  is  predominant  wave  direction.  This  certainly  would 
be  true  for  sites  naturally  protected  in  all  but  one  narrow  sector.  WIS 
Phase  III  data  did  not  include  any  such  sites,  however,  assuming  an  open  coast 
with  sheltering  only  from  major  landforms. 

Figure  7  illustrates  the  time  series  for  wave  heights  during  October 
1956  at  Nagshead,  North  Carolina.  This  particular  time  span  was  chosen  for 
presentation  because  it  Included  rapid  changes  in  wave  conditions,  especially 
on  October  27  and  28,  1956,  as  indicated  by  sharp  spikes  near  the  end  of  the 
wave  height  time  series  plot  of  Figure  7.  Table  Al  includes  Phase  III  wave 
information  recorded  for  these  2  days. 

Figure  A6  shows  the  time  series  of  peak  wave  period  during  this  same 


Figure  7.  Wave  height  time  series:  Nagshead,  North  Carolina, 

October  1956 


month  for  Nagshead,  North  Carolina.  The  wave  period  can  be  seen  to  vary  some¬ 
what  out  of  phase  with  wave  height  and  to  have  a  tendency  to  remain  constant 

2 

for  significant  time  spans  and  then  change  abruptly.  A  plot  of  H/gT  for 

the  same  period  at  Nagshead  (Figure  A7)  appears  more  like  the  wave  period  time 

series  than  the  wave  height  time  series,  also  tending  to  vary  slowly  for  sig- 

2 

nif leant  time  spans  and  change  abruptly  (the  influence  of  T  in  the  denomi- 

2 

nator).  A  plot  of  wave  severity,  H  L  (Figure  A8) ,  for  the  same  time  period 

dramatically  delineates  extremes  of  the  wave  height  time  series.  When  plotted 

2  ]  /  3 

in  Figure  A9  as  (H  L)  ,  however,  wave  severity  very  closely  resembles  the 

wave  height  time  series  plot.  Wave  severity  in  this  form  has  the  same  units 

2 

as  wave  height  and  includes  the  influence  of  H  to  balance  the  influence  of 
in  the  denominator.  The  plot  of  direction  (Figure  A10)  does  not  indicate 
direct  relation  to  the  wave  height  plot  and  is  much  more  erratic,  even  in  non¬ 
extreme  periods.  Direction  can  be  considered  to  be  practically  independent  of 
the  intensity  of  wave  conditions  since  it  is  controlled  almost  exclusively  by 
geometric  factors. 

The  convention  of  previous  investigators  (Houmb  and  Vik  1977)  to  rely 
solely  on  variation  in  zero  moment  wave  height  for  definition  of  extreme  event 

i 

i 
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durations  was  maintained  in  this  study.  This  parameter  is  most  easy  to  visu¬ 
alize  and  has  a  long  tradition  as  the  critical  measure  of  Intensity  of  extreme 

2  2  1/3 

events  at  sea.  Variations  of  H  L  and  (H  L)  show  promise,  but  the  large 
2 

units  of  H  L  make  results  of  computations  rather  abstract  and  the  variation 
2  1/3 

of  (H  L)  seemed  quite  close  to  that  of  H  .  Relationships  of  individual 
extreme  event  durations  (measured  by  variation  of  H)  to  peak  conditions  mea¬ 
sured  by  all  parameters  discussed  above  were  investigated,  however,  and  the 
results  of  that  analysis  are  reported  later  in  this  report. 

An  investigation  of  actual  weather  conditions  on  the  Atlantic  coast  in 
the  time  frame  surrounding  October  27-28,  1956,  was  conducted  to  better  under¬ 
stand  what  events  were  actually  driving  the  numerical  simulations  to  produce 
irregularities  in  the  time  series  of  Figure  7.  First,  Phase  II  data  input  to 
the  Phase  III  numerical  wave  transformation  were  inspected.  Table  A2  presents 
Phase  li  information  at  Station  A2037,  at  36.06°  N  latitude  and  74.92°  W  lon¬ 
gitude,  approximately  33  nautical  miles  (61  km)  east-southeast  of  Nagshead  in 
about  240  ft  (73  ml  ot  water.  The  intermittent  appearance  and  disappearance 
of  swell  can  be  seen  to  follow  a  similar  pattern  in  the  Phase  III  site  of 
interest  (Station  A3083)  and  the  Phase  II  site  directly  offshore  (Sta¬ 
tion  A2037).  Wave  heights  In  deeper  water  are  higher,  lacking  the  depth  limi¬ 
tations  inherent  in  Phase  HI  simulations.  Wave  periods  of  both  sites  are 
identical,  unaffected  by  the  wave  transformation  processes  simulated  in 
Phase  III.  The  direction  convention  in  Phase  II  is  different,  indicating  the 
direction  from  which  waves  are  travelling  toward  the  center  of  the  compass 
rose.  Phase  II  data  do  not  include  anything  significantly  revealing  about  the 
irregularities  of  interest,  basically  showing  the  same  patterns  in  this  case. 
The  nearest  Phase  I  site  offshore  of  Nagshead  was  Station  A1005  at 


35.4°  N  latitude  and  72.3°  W  longitude,  located  in  deep  water  approximately 
163  nautical  miles  (302  km)  east-southeast  of  Nagshead,  North  Carolina. 

Table  A3  shows  Phase  I  information  recorded  for  October  26-28,  1956.  There  is 
only  one  record  which  included  swell;  that  record  did  not  dominate  the  com¬ 
bined  wave  height,  which  appears  to  be  steadily  decreasing  at  that  time.  It 
is  Important  to  recognize  that  a  significant  travel  time  would  be  involved 
between  this  Phase  I  site  and  Stations  A2037  or  A3083  (approximately  8  and 
10  hr,  respectively,  for  waves  of  11-sec  period),  so  the  conditions  at  a  given 
date  and  time  should  be  "out  of  phase"  by  three  to  four  records. 

An  inspection  of  surface  weather  charts  during  the  later  part  of  October 
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1956  for  North  America  and  the  north  Atlantic  Ocean  was  made  to  identify  syn¬ 
optic  weather  systems  which  may  have  dominated  Phase  I  and  Phase  II  informa¬ 
tion.  An  explanation  was  sought  for  the  sudden  appearance  and  disappearance 
of  swell  in  the  data,  as  well  as  an  explanation  of  differences  between  Phase  I 
and  Phase  II  wave  information.  Figure  8  illustrates  recorded  weather  patterns 
of  October  26-28,  1956,  showing  the  presence  of  a  generally  stationary,  weakly 
defined,  low  pressure  system  of  fluctuating  intensity  offshore  of  Cape 
Hatteras.  This  location  is  close  to  Station  A1005;  thus,  the  basic  definition 
of  swell  as  waves  which  have  left  their  area  of  generation  could  explain  the 
lack  of  swell  in  Phase  I  data.  The  wave  field  at  this  point  would  have  been 
under  the  influence  of  cyclonic  winds  of  the  low  pressure  system  and  thus  only 
seas  would  have  existed,  as  defined  by  WIS  conventions.  The  relative  position 
of  Stations  A2037  and  A3083  in  combination  with  the  fluctuating  intensity  of 
the  low  pressure  system  appears  to  have  caused  swell  either  to  come  from  too 
tar  south  to  affect  Nagshead  or  to  exist  only  as  seas,  except  for  the  spikes 
of  Figure  7.  This  set  of  circumstances  is  probably  exceptional,  but  an  under¬ 
standing  of  the  real  weather  patterns  driving  numerical  simulations  of  the  WIS 
program  in  this  instance  may  help  explain  trends  of  duration  revealed  by  fur¬ 
ther  analysis  of  WIS  data. 


Analytical  Procedure  and  Results 


A  FORTRAN  computer  program  was  written  which  read  the  58,440  records 
stored  for  each  Phase  III  site  and  maintained  a  record  of  the  number  of  consec¬ 
utive  records,  each  of  which  had  a  combined  wave  height  above  a  specified  wave 
height  threshold,  HI  .  Subsequent  us*  of  the  term  "extreme  event"  refers  to 
events  defined  in  this  manner.  The  number  of  extreme  events  was  counted  and 
statistics  including  the  maximum,  minimum,  mean,  and  standard  deviation  dura¬ 
tions  were  computed.  Peak  conditions  of  each  extreme  event  were  noted  as  the 
highest  combined  wave  height  in  a  consecutive  series  above  the  threshold,  and 
the  period  and  direction  of  sea  or  swell,  whichever  had  the  highest  incremen¬ 
tal  wave  height.  Maximum,  minimim,  mean  and  standard  deviation  wave  heights 
also  were  computed.  Each  data  set  included  20  years  of  record,  so  the  number 
of  extreme  events  per  year  (the  Poisson  lambda  parameter!  was  computed  as  the 
total  number  of  events  divided  by  20. 

Initial  runs  of  this  extreme  event  identification  program  resulted  in  a 
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Figure  8.  Eastern  US  surface  weather  patterns:  October  1956 


surprisingly  large  number  of  extreme  events,  consistently  on  the  order  of 
30  to  40  percent  of  all  extreme  events  identified,  to  be  only  3  hr  in  dura¬ 
tion,  i.e.,  only  one  record  above  the  threshold.  The  actual  duration  could  be 
anywhere  from  0  to  6  hr  for  a  single  record  above  the  threshold,  but  an  aver¬ 
age  value  of  3  hr  was  consistently  assumed  in  such  cases.  Variation  of  the 
threshold  had  little  effect  on  the  percentage,  although  the  total  number  of 
extreme  events  was  of  course  affected.  Review  of  climatology  considerations 
inherent  in  WIS  simulations  (Kesio  and  Hayden  1973;  and  Corson,  Reslo,  and 
Vincent  1980)  did  not  uncover  a  rationale  for  excluding  a  priori  durations 
that  short.  In  fact,  a  duration  of  3  hr  is  either  implicitly  or  explicitly 
assumed  for  peak  conditions  in  many  wave  forecasts,  designs,  and  research 
efforts  relating  to  the  tradition  of  sampling  wave  gages  at  this  interval. 
Average  low  pressure  systems  which  would  generate  extreme  wave  conditions  are 
known  to  typically  last  much  longer,  however,  as  in  the  case  of  the  system 
illustrated  in  Figure  8. 

In  view  of  this  last  fact  and  the  example  of  late  October  1956  at 
Nagshead,  the  program  was  adjusted  to  ignore  a  lapse  below  the  threshold  of 
only  one  record  (i.e.,  6  hr)  between  consecutive  extreme  events,  as  Identified 
previously.  This  adjustment  lowered  the  number  of  extreme  events  of  only  3  hr 
duration  (one  record  above  the  threshold)  only  slightly,  but  a  neglect  of 
longer  lapses  or  other  adjustments  to  the  identification  procedure  could  not 
be  rationalized.  Tables  A4-A8  give  duration  results,  following  the  procedures 
described  above,  for  the  five  sites  at  all  thresholds  investigated. 

Mean  and  maximum  durations  for  the  three  Atlantic  sites  are  virtually 
identical  to  those  reported  by  Jensen  (1983a),  with  the  occasional  exception 
caused  by  combination  of  two  events  separated  by  only  one  record  with  H 
below  the  threshold.  The  mean  duration  was  slightly  higher  in  these  few 
cases . 


The  percent  occurrence  of  wave  heights  (percent  records  H  >  HI)  was  of 
special  interest  since  this  statistic  for  a  range  of  HI  levels  is  now  or 
will  be  published  and  readily  available  for  all  WIS  stations  of  all  three 
phases.  It  was  hoped  this  nondimensional  parameter  could  be  used  as  a  tool 
for  choosing  threshold  levels  for  duration  computations  which  would  preclude 
manv  of  the  iterations  which  otherwise  might  be  necessary.  The  number  of  ex¬ 
treme  events  per  year  was  also  of  special  interest  since  this  parameter  is  so 
important  in  extremal  statistics  and  expectations. 
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Figure  9  shows  the  mean  and  standard  deviation  of  duration  plotted 
against  percent  occurrence  (actually  exceedance)  of  wave  heights  above  the 


Nagshead.  North  Carolina 


Figure  9.  Mean  duration  and  standard  deviation  versus  percent 
occurrence  of  wave  height  threshold,  Nagshead,  North  Carolina 


specified  threshold  for  Nagshead.  Figures  A11-A14  show  mean  and  standard  de¬ 
viation  durations  for  the  other  four  sites  plotted  against  percent  occurrence 
of  wave  heights  above  the  threshold.  Figure  10  shows  the  nearly  linear  rela¬ 
tionship  of  the  number  of  extreme  events  per  year  with  percent  occurrence  of 
wave  heights  above  a  specified  threshold  for  Nagshead.  A  similar  trend  is 
evident  for  higher  wave  with  similar  percent  occurrence  at  Newport,  Oregon,  as 
shown  in  Figure  A15.  These  plots  in  themselves  do  not  indicate  an  outstanding 
range  of  percent  occurrence  as  a  choice  for  definition  of  extreme  events  and 
durations.  Some  subjective  choices  can  be  made  since  an  important  purpose  of 
this  exercise  is  to  identify  extreme  events.  Clearly,  an  excessively  large 
number  of  extreme  events  per  year,  say  more  than  20,  will  probably  include 
some  events  that  can  hardly  be  regarded  as  "extremes"  in  the  practical  sense. 
On  the  other  hand,  an  average  number  of  extreme  events  per  year  less  than  one 
or  two  would  generally  imply  exclusion  of  some  events  which  belong  in  a 
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Figure  10.  Extreme  events  per  year  versus  percent  occurrence, 

Nagshead,  North  Carolina 

population  of  extremes.  These  considerations  are  consistent  with  the  author  s 
experience  in  developing  design  criteria  based  on  extremal  statistics  of  peak 
wave  height  conditions  (e.g.  Andrew,  Smith,  and  McKee  19851. 

A  simple  linear  regression  of  extreme  events  per  year  with  percent  oc¬ 
currence  of  wave  heights  above  the  threshold,  constrained  to  pass  through  the 
origin,  for  the  41  cases  considered  at  all  five  sites  indicates  that  percent 
occurrence  *  0.3X  with  a  correlation  coefficient  of  0.97.  This  relation  ap¬ 
plies  to  both  the  Atlantic  and  Pacific  sites  addressed  individually,  even 
though  the  absolute  value  of  wave  heights  themselves  on  the  Pacific  are  sub¬ 
stantially  higher  than  those  on  the  Atlantic  at  the  same  percent  occurrence 
levels.  A  range  in  X  of  2  to  20  thus  would  correspond  to  a  range  in  percent 
occurrence  of  0.6  to  6.0  percent  for  the  choice  of  a  desirable  threshold 
level,  HI  .  The  lower  limit  of  this  range  would  guarantee  a  sample  size  of 
at  least  40  extreme  events,  which  is  generally  desirable  for  most  statistical 
considerations.  The  choice  of  a  threshold  wave  height  may  be  made  more  pre¬ 
cisely  when  some  physical  tolerance  level  is  at  issue,  for  example  the  point 
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at  which  some  operation  at  sea  must  be  temporarily  terminated. 

The  other  parameters  presented  in  Tables  A4-A8  show  interesting  trends. 
The  minimum  duration  was  3  hr  in  every  case  except  one  where  only  three  ex¬ 
treme  events  were  identified.  A  count  of  extreme  events  with  a  3-hr  duration 
for  the  Nagshead  cases  indicated  32  to  48  percent  of  extreme  events  shared  the 
minimum  duration.  No  relation  of  the  number  of  extreme  events  with  a  3-hr  du¬ 
ration  to  the  threshold  level  was  apparent.  The  maximum  duration  can  be  seen 
to  be  proportional  to  percent  occurrence  of  wave  heights  above  the  threshold 
and  typically  many  standard  deviations  above  the  mean.  The  mean  duration  ac¬ 
cordingly  also  is  proportional  to  percent  occurrence  of  wave  heights  above  the 
threshold.  The  standard  deviation  was  rarely  less  than  the  mean,  but  always 
of  the  same  order  of  magnitude.  A  lack  of  central  tendency  for  durations  was 
noted  by  Houmb  and  Vik  (1977). 

Another  scheme  of  extreme  event  identification  was  investigated  which 
actually  applied  a  lower  threshold  HI  in  the  same  way  for  determination  of 
duration,  but  only  to  extreme  events  whose  peak  (combined)  wave  height  was 
above  a  second  higher  threshold,  H2  .  The  most  notable  effects  of  the  second 
threshold  were  to  substantially  reduce  the  number  of  extreme  events  per  year 
for  a  given  HI  threshold  and  to  reduce  the  number  of  extreme  events  with  a 
3-hr  duration  to  zero  in  nearly  every  case.  Variation  of  HI  with  a  fixed 
H2  had  little  effect  on  the  number  of  extreme  events  per  year.  The  central 
tendency  of  durations  was  somewhat  stronger  in  these  subsets,  with  the  stan¬ 
dard  deviation  often,  but  not  always,  less  than  the  mean.  These  two  parame¬ 
ters  consistently  retained  the  same  order  of  magnitude. 

Tables  A9  and  A10  present  the  parameter  values  computed  for  various 
combinations  of  HI  and  H2  at  the  Nagshead  and  Daytona  Beach  sites.  Fig¬ 
ures  A16  and  A17  show  variation  of  the  mean  and  standard  deviation  durations 
with  percent  occurrence  of  the  lower  threshold  HI  at  a  upper  threshold  H2 
tor  peak  conditions  fixed  at  300  cm  (0.6  percent)  and  300  cm  (0.8  percent)  for 
the  Newport  and  Nagshead  sites.  This  scheme  of  double  thresholds  for  extreme 
event  identification  was  not  pursued  further  since  it  was  considered  more  de¬ 
sirable  to  address  trends  in  peak  conditions  separately  from  durations  above  a 
specified  threshold.  An  approach  which  addressed  marginal  distributions  ver¬ 
sus  conditional  distributions  was  preferred. 


CHAPTER  VII:  DISTRIBUTION  OF  DURATIONS 


Method  of  Analysis 

The  cumulative  probability  of  durations  derived  by  the  single  threshold 
method  described  above  is  estimated  by  application  of  a  plotting  formula  corn- 
monly  applied  in  analyses  of  this  type  (Gumbel  1958,  and  Isaacson  and 
MacKensie  1981): 


F(tt)  =  ~  Y)  *  =  1.  2,  3 . n  (69) 

where  F(t^)  Is  the  estimated  cumulative  probability  of  the  "it^"  smallest 
duration  and  n  is  the  number  of  extreme  events.  Durations  are  first  ordered 
from  smallest  to  largest  for  this  purpose  and  the  corresponding  cumulative 
probability  computed.  Other  plotting  formulae  were  considered  (e.g., 
Gringorten  1963)  ,  but  this  more  commonly  used  approach  is  preferable  for  gen¬ 
eral  application  since  no  additional  parameters  need  be  estimated. 

Two  continuous  distributions  are  considered  as  models  for  the  cumulative 
probability  of  durations  because  of  their  common  application  to  peak  wave 
height  conditions:  the  Extremal  (Fisher-Tippett)  Type  I  and  the  Weibull  dis¬ 
tributions.  An  existing  FORTRAN  program  (US  Army  Engineer  Waterways  Experi¬ 
ment  Station  1985)  ,  originally  designed  to  fit  these  distributions  to  wave 
height  data  by  the  method  of  least  squares,  applying  the  plotting  formula  con¬ 
vention  described  above,  was  adapted  to  work  instead  with  duration  data 
derived  by  the  extreme  event  identification  program.  The  extreme  event  iden¬ 
tification  program  was  ultimately  combined  with  the  program-estimating  distri¬ 
bution  parameters  and  titled  STRMDIST,  a  listing  of  which  is  presented  in 
Appendix  C.  The  program  STRMDIST,  in  addition  to  the  extreme  event  indentifi- 
cation  and  duration  derivation  computations  already  described,  computes  dis¬ 
tribution  parameters  (e  and  8  for  the  Extremal  Type  I  and  a  and  B  for 
the  Weibull),  estimated  (distribution)  mean  and  standard  deviation,  correla¬ 
tion  coefficient,  sum  of  the  square  residuals,  and  standard  error.  These  pa¬ 
rameters  also  are  computed  for  peak  wave  heights  of  extreme  events  identified. 
Tables  A 1 1 —A  1 5  give  results  of  the  STRMDIST  analysis  for  five  Phase  III  sites. 
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Discussion  of  the  Distribution  Analysis 
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Figures  11  and  12  demonstrate  fit  of  the  least  squares  regression  dis¬ 
tribution  to  the  data  as  represented  by  the  plotting  formula  for  one  case  each 

Nagshead.  North  Carolina  -  Threshold  H  -  300  cm 


Cumulative  Probability 

Figure  11.  Duration  cumulative  probability: 

Nagshead,  North  Carolina 

at  Nagshead,  North  Carolina,  and  Newport,  Oregon.  The  Weibull  distribution  in 
both  these  cases  can  be  seen  to  generally  fit  the  overall  data  spread  better, 
but  the  Extremal  Type  I  comes  closer  to  the  few  most  extreme  durations.  The 
correlation  coefficients,  sums  of  square  residuals,  and  standard  errors  in 
Tables  All —A 1 5  indicate  that  the  Weibull  distribution  generally  fits  the  data 
better  than  the  Extremal  Type  I,  but  both  distributions  fit  it  acceptably  well 
in  practical  terms.  Correlation  coefficients  above  0.90  would  provide  a  rule- 
of-thumb  acceptable  fit  in  exercises  of  this  type  with  weather-related  data. 
Both  distributions  generally  exceed  this  criterion. 

Figure  13  shows  correlation  coefficients  for  both  distributions  plotted 
against  percent  occurrence  of  wave  heights  above  the  specified  threshold, 

HI  ,  for  Nagshead.  Figure  A16  shows  a  similar  plot  for  Newport. 
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Figure  12.  Duration  cumulative  probability:  Newport,  Oregon 
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Nagshsad,  North  Carolina 


Percent  Occurrence  of  Threshold  Wave  Height 


Figure  13.  Correlation  coefficient  versus  percent  occurrence 
of  wave  height  threshold,  Nagshead,  North  Carolina 


Figures  A17  and  A18  show  the  correlation  coefficients  for  both  distributions 
at  Nagshead  and  Newport  plotted  against  the  number  of  extreme  events  per  year. 
No  obvious  maximum  occurs  which  could  reliably  be  taken  as  an  indication  of  an 
optimal  choice  for  either  X  or  HI  . 

Figures  A19  and  A21  are  plots  of  the  standard  error  against  extreme 
events  per  year  and  percent  occurrence  for  Nagshead,  North  Carolina.  Fig¬ 
ures  A20  and  A22  show  the  same  information  for  Newport,  Oregon.  Again,  no  ob¬ 
vious  minimum  generally  occurs  to  indicate  an  optimal  choice  for  X  or  HI  . 

Figures  14  and  15  are  graphs  of  the  sample  and  distribution  means  and 
sample  and  distribution  standard  deviations  plotted  against  percent  occurrence 
and  the  number  of  extreme  events  per  year  both  for  Nagshead.  Figures  A23 
and  A24  are  similar  graphs  for  Newport.  The  Extremal  Type  I  distribution  mean 
and  standard  deviation  can  be  seen  to  generally  come  closer  to  the  sample  mean 
and  standard  deviation.  This  is  desirable,  particularly  in  the  case  of  the 
mean.  The  Central  Limit  Theorem  states  that  sample  means  from  an  infinite 
population  can  be  considered  as  random  variables  with  a  mean  equal  to  the  pop¬ 
ulation  mean.  The  standard  deviation,  as  a  measure  of  the  spread  of  duration 
values  about  the  mean,  is  an  important  indicator  of  how  conservative  a  param¬ 
eterized  distribution  might  be.  The  Extremal  Type  I  distribution  can  be  seen 
to  be  closer  to  and  consistently  larger  than  (i.e.  on  the  conservative  side 
of)  the  sample  standard  deviation.  The  Weibull  distribution  standard  devia¬ 
tion  is  both  farther  from  the  sample  standard  deviation  and  generally  lower, 
i.e.,  predicting  more  central  tendency  than  the  sample.  The  Extremal  Type  I 
distribution  in  these  respects  appears  superior  to  the  Weibull  distribution. 


50 


Duration  Standard  Deviation  (hours) 


Nagahead.  North  Carolina 


Figure  14.  Mean  duration  versus  percent  occurrence  of 
wave  height  threshold,  Nagshead,  North  Carolina 
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Percent  Occurrence  of  Threshold  Wave  Height 

Figure  15.  Duration  standard  deviation  versus  percent 
occurrence  of  wave  height  threshold,  Nagshead, 

North  Carolina 


CHAPTER  VIII:  RELATIONSHIP  OF  DURATION  TO  PEAK  CONDITIONS 


The  potential  linear  or  nonlinear  relationship  of  an  extreme  event's 
duration  with  peak  conditions  of  the  extreme  event  were  investigated  with  the 
aid  of  statistical  software  package  SPSS  (Nie  et  al.  1975).  The  stepwise  mul¬ 
tiple  regression  capabilities  of  SPSS  were  of  particular  value  in  testing 

whether  extreme  event  duration  appeared  to  be  dependent  on  peak  conditions,  as 

2  2  2  2 

measured  by  various  parameters  such  as  H,  T,  H  ,  T  ,  H/gT  ,  HL, 
and  direction.  Simple  linear  regressions  of  extreme  event  durations,  as  de¬ 
rived  by  a  range  of  thresholds,  first  were  performed.  In  the  same  program  ex¬ 
ecution,  SPSS  allowed  a  stepwise  multiple  regression  of  duration  against  H  , 

2  2 

H  ,  T  ,  and  T  to  be  performed.  This  procedure  estimated  the  incremental 
contribution  of  each  of  these  potentially  controlling  (independent)  variables 
to  the  data  fit  by  the  least  squares  method.  An  equation  of  the  following 
form  was  thus  possible,  assuming  the  contribution  of  each  of  these  tested 
parameters  was  significant: 

t  =  aH  +  bH2  +  cT  +  dT2  (70) 

where  a  ,  b  ,  c  ,  and  d  are  constants. 

The  purpose  in  this  exercise  was  not  to  derive  a  predictor  equation,  but 

to  see  if  a  significant  relationship  existed.  Therefore,  the  obvious  inter- 
2  2 

dependence  of  H  with  H  and  T  with  T  was  not  of  undesirable  conse¬ 
quence.  One  common  technique  to  test  for  existence  of  a  nonlinear  relation¬ 
ship,  versus  a  linear  relationship,  is  also  to  test  the  square  of  the  variable 
on  a  trial  basis.  A  substantially  improved  fit  with  the  square  of  the  parame¬ 
ter  included  in  the  regression  equation  generally  indicates  that  a  nonlinear 
relationship,  whether  polynomial  or  otherwise,  is  more  reliable  than  a  simple 
linear  relationship.  The  correlation  coefficient,  r  ,  as  applied  above  in 
the  fit  of  distribution  functions,  was  taken  as  the  primary  measure  of  the 
strength  of  a  relationship  in  this  analysis. 

Tables  A16-A20  show  results  from  execution  of  SPSS  for  all  cases  tested 
for  each  of  the  five  Phase  III  sites.  A  listing  of  the  SPSS  command  file  used 
to  perform  each  of  these  executions  is  presented  in  Appendix  D  along  with  a 


sample  output.  The  tables  give  correlation  coefficients  for  duration  against 
2  2 

H  ,  H  ,  T  ,  and  T  (individual  simple  linear  regressions) ,  against  all 

2 

four  of  these  parameters  in  a  stepwise  procedure  and  against  H  L  (simple 

linear  regression).  There  is  little  indication  in  any  case  of  a  linear  rela- 

2 

tionship  of  duration  with  H/gT  with  correlation  coefficients  for  this  pa¬ 
rameter  consistently  near  zero.  Similarly,  correlation  coefficients  of  dura¬ 
tion  with  predominant  direction  of  wave  propagation  at  the  peaks  of  extreme 
events  were  consistently  near  zero. 


Discussion  of  Results  of  the  Regression  Analysis 


The  parameter  h  ,  the  peak  zero  moment  wave  height,  is  consistently  the 
most  significant  parameter,  which  confirms  that  an  extreme  event  identifica¬ 
tion  and  duration  definition  procedure  using  this  parameter  is  best.  Another 
notable  trend  indicated  by  the  above  results  is  the  observation  that  correla¬ 
tion  coefficients  for  the  Pacific  sites  are  consistently  lower  than  those  for 
Atlantic  sites.  A  possible  explanation  for  this  is  the  fact  that  Pacific 
storms  typically  form  well  away  from  the  coast  and  travel  onshore.  They  tend 
to  be  well  formed  when  their  effects  first  become  significant  and  their  tracks 
are  more  or  less  in  the  same  direction  (eastward  to  some  degree).  Atlantic 
(extratropical)  storms  can  form  onshore  and  travel  seaward,  travel  longshore, 
or  linger  in  one  spot,  as  exemplified  by  the  previous  account  of  conditions  in 
late  October  1956.  This  more  variable  track  (particularly  the  potential  for  a 
roughly  stationary  storm)  may  cause  the  duration  above  a  specified  threshold 
in  many  cases  to  be  more  dependent  on  the  time-history  of  the  storm’s  internal 
intensity  than  its  track  past  a  fixed  site. 

There  was  no  strong  correlation  of  duration  (applying  the  rule-of-thumb 
criterion  of  0.90)  with  any  of  the  variables  on  either  coast.  The  regression 
slopes,  i.e.  the  B  parameter  in  Equation  4U,  also  were  consistently  small 
numbers,  much  closer  to  zero  than  to  one.  The  low  slopes,  even  for  H  ,  indi¬ 
cate  that  dependence  of  durations  on  peak  wave  conditions  is  weak.  A  fully 
rigorous  proof  of  dependence  or  otherwise  would  require  many  more  tests  and 
computations  than  those  presented  here.  The  lack  of  an  obvious  strong  depen¬ 
dence,  however,  raises  the  suggestion  that,  for  practical  purposes,  extreme 
event  duration  might  be  taken  as  independent  of  peak  conditions  of  the  extreme 
event.  This  would  make  estimates  of  joint  probability,  for  example  forecast 


of  durations  of  wave  heights  above  a  threshold  for  a  rare  event  (e.g.,  the  50- 
or  100-year  extreme  event),  relatively  easy  to  compute.  An  example  of  how 
such  an  estimate  might  be  made  follows: 


Example  Computation  of  Peak  Wave  Height  and  Duration 
Joint  Probability 


Problem:  What  is  the  joint  probability  of  zero  moment  wave 
heights  greater  than  3.0  m  lasting  longer  than  12  hr  during  an 
extreme  event  whose  peak  zero  moment  wave  height  is  greater  than 
4.5  m  at  Nagshead,  North  Carolina? 

Solution:  The  definition  of  duration  at  HI  allows  the  asso¬ 

ciated  parameters  presented  in  Table  A12  to  be  applied.  Choosing 
the  Extremal  Type  I  distribution  to  represent  both  marginal  dis¬ 
tribution  of  peak  wave  heights  and  marginal  distribution  of  dura¬ 
tions:  e  =6.30,  $  =15.8,  e„  =  326.3  ,  and  R  =  48.0  . 

t  t  H  H 

The  Poisson  parameter,  X  ,  from  Table  A5,  is  3.8.  The  marginal 
probabilities  of  exceedance  are: 


=  0.073 


The  joint  probability,  taken  as  the  product  of  independent  mar¬ 
ginal  probabilities  defined  from  the  same  population  (HI  =  300), 
is : 


P(t’  >  12  ,  H'  >  450  |  HI  -  300)  =  0.502(0.073)  =  0.037 


The  associated  return  period  Is: 


RT(t>H)  ■  rm  -  F'a.Hm  ■  ■  7-° yea” 


The  associated  nonencounter  probability  in  a  50-year  time  period 
is : 


NE(t,H)  =  exp  -  exp  *  0.00079  -  0.08% 


The  associated  risk  of  encountering  such  a  condition  in  a  50-year 
time  span  is  1  -  NE(t,H)  »  0.921  »  92.1%. 

Discussion:  Given  the  assumptions  stated  above,  the  probability 
of  exceedance  of  a  peak  wave  height  of  4.5  m  of  any  duration  is 
7.3  percent.  The  condition  of  duration  exceeding  12  hr  eliminates 
about  half  of  the  possibilities;  therefore,  the  joint  probability 
is  about  half  as  much.  The  joint  return  period  is  also  corre¬ 
spondingly  longer.  The  Poisson  assumption  inherent  in  definition 
of  return  period  and  nonencounter  probability  can  be  extended  to 
the  joint  peak  wave  height  and  duration  distribution  if  waiting 
periods  between  extreme  events  are  much  greater  than  durations  of 
the  extreme  events.  The  Poisson  distribution  is  a  discrete  dis¬ 
tribution,  and  its  application  technically  extends  only  to  dis¬ 
crete  events. 
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CHAPTER  IX:  CONCLUSIONS 


Literature  Review 

A  review  of  scientific  and  engineering  literature  related  to  duration  of 
sea  states  reveals  little  direct  work  in  this  area.  The  work  of  Houmb  and  Vik 
(1977)  is  most  pertinent  to  objectives  of  this  study.  These  investigators 
worked  with  several  years  of  intermittently  measured  wave  information  at  five 
points  along  the  North  Sea  coast  of  Norway.  They  found  the  duration  of  ex¬ 
treme  sea  states,  as  defined  by  the  exceedance  of  a  wave  height  threshold,  to 
fit  a  Weibull  distribution.  They  approached  the  problem  as  much  as  possible 
from  a  theoretical  perspective  in  order  to  maximize  the  reliability  of  obser¬ 
vations  based  on  limited  data. 

Identification  of  Extreme  Events 


This  study  applies  the  Phase  III  (shallow  water)  Wave  Information  Stud¬ 
ies  (WIS)  database  of  hindcast  wave  data  because  of  its  unusually  long,  con¬ 
tinuous  20-year  period  of  record  and  because  of  its  synoptic  (ocean  wide) 
perspective  on  wave  conditions.  The  WIS  numerical  simulations  involve  some 
practical  simplifications,  but  no  database  of  measured  wave  information  is 
available  which  could  be  used  to  investigate  such  a  long  period  of  record  over 
a  wide  geographical  area.  Data  from  five  Phase  III  stations  are  applied  in 
this  study  to  investigate  duration  of  extreme  wave  conditions.  Three  are  on 
the  Atlantic  coast  (from  New  Jersey  to  central  Florida)  and  two  are  on  the 
Pacific  coast  (Oregon  to  central  California). 

The  conventional  parameter  for  long-term  wave  statistics,  zero  moment 
wave  height,  is  chosen  as  the  most  practical  and  reliable  indicator  of  inten¬ 
sity  of  wave  conditions.  A  computer  program  is  presented  which  reviews 
Phase  III  information  and  records  the  number  of  sequential  records  (each  3  hr 
apart)  in  which  the  geometric  average  (combined)  sea  and  swell  wave  height  is 
above  a  specific  threshold.  A  single  record  below  the  threshold  between  two 
that  were  above  is  ignored,  i.e.,  the  two  records  above  are  treated  as  part  of 
a  single  event.  The  percent  occurrence  of  waves  above  a  threshold  is  found  to 
vary  linearly  with  the  number  of  extreme  events  identified,  regardless  of  ab¬ 
solute  intensity  of  wave  climate  on  either  coast. 


Distribution  of  Durations 


The  Weibull  and  Extremal  Type  I  distributions  are  fit  by  the  method  of 
least  squares  to  durations  of  extreme  events  identified  and  to  peak  wave 
heights.  Both  distributions  show  acceptable  correlation  to  the  wave  data,  but 
the  Extremal  Type  I  is  found  to  provide  superior  estimates  of  both  durations 
and  peak  wave  heights. 

Relationship  of  Duration  to  Peak  Intensity 

A  multilinear  regression  analysis  is  performed  to  address  the  potential 
relationship  of  extreme  event  duration  to  peak  conditions  of  the  extreme 
event.  Peak  intensity,  as  measured  by  the  zero  moment  wave  height,  has  only  a 
weak  linear  relationship  to  duration.  Other  alternate  parameters  of  intensity 
show  little  evidence  of  significant  linear  relation  to  duration.  The  investi¬ 
gation  does  not  rigorously  prove  statistical  independence,  but  the  assumption 
of  independence  of  duration  from  peak  intensity  is  proposed  as  an  expedient 
measure.  This  assumption  greatly  simplifies  prediction  of  durations  of  wave 
conditions  above  a  critical  threshold. 
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Figure  A2.  WIS  Phase  I  grid.  North  Pacific  Ocean 
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Figure  A4.  WIS  Phase  II  grid.  Pacific  coast 
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Figure  A5.  Mid-Atlantic  coast  portion,  WIS  Phase  III  stations 
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Figure  A6.  Wave  period  time  series,  Nagshead,  NC,  October  1956 
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Figure  A7 .  Wave  steepness  time  series,  Nagshead,  NC,  October  1956 
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rigure  A8.  Wave  severity  time  series,  Nagshead,  NC,  October  1956 
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Figure  A9.  Wave  severity  time  series,  Nagshead,  NC,  October  1956 
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Figure  A10.  Wave  direction  time  series,  Nagshead,  NC,  October  1956 
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Figure  A12.  t  and  st  versus  percent  occurrence, 

Newport,  OR 
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Figure  A15.  t  and  st  versus  percent  occurrence  with 
peak  threshold,  Daytona  Beach,  FL 
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Figure  A21.  Standard  error  versus  extreme  events  per  year, 

Nagshead,  NC 
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Figure  A22.  Standard  error  versus  extreme  events  per  year, 

Newport,  OR 
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Figure  A23.  Mean  duration  versus  percent  occurrence, 

Newport,  OR 
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Figure  A24.  Duration  standard  deviation  versus 
percent  occurrence,  Newport,  OR 


Table  A1 

October  1956  Phase  III  Data,  Nagshead,  NC 
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Table  A4 


Number  Number 

HI  H  >  HI  H  >  HI  of  of  tmin  fcmax  t  °t 

cm  cm  %  Events  Events/yr  hrs  hrs  hrs  hrs 


200  1442 

250  384 

300  81 

350  18 


2.5 

323 

16.2 

3 

54 

12.8 

10.9 

0.7 

112 

5.6 

3 

30 

9.9 

7.3 

0. 1 

29 

1.4 

3 

24 

8.1 

6.8 

0.03 

9 

0.4 

3 

15 

5.7 

4.1 

Table  A5 

Duration  Information  for  Nagshead,  NC 


A22 


HI 

cm 

H  >  HI 

cm 

H  >  HI 

t 

Number 

of 

Events 

Number 

of 

Events/yr 

fcmin 

hrs 

tmax 

hrs 

t 

hrs 

°t 

hrs 

150 

6983 

11.9 

79  2 

39.6 

3 

570 

26.9 

39.1 

200 

3167 

5.4 

460 

23.0 

3 

306 

21.1 

29.4 

250 

1093 

1.9 

179 

9.0 

3 

165 

18.9 

25.0 

300 

374 

0.6 

77 

3.8 

3 

111 

15.1 

17.2 

350 

143 

0.2 

36 

1.8 

3 

84 

12.2 

14.9 

400 

56 

0. 10 

13 

0.6 

3 

42 

12.9 

12.0 

450 

16 

0.03 

8 

0.4 

3 

15 

6.4 

4.1 

I 

i. 

! 


i 

i 


Table  A6 

Duration  Information  for  Davtona  Beach,  FL 


5 

i  hi 

H  >  HI 

H  >  HI 

Number 

of 

Number 

of 

u  S 

cm 

?  .. 

Events 

Events/vr 

150 

8855 

1' 

200 

4183 

i  250 

< 

1340 

j  300 

478 

i 

-  350 

r 

143 

i 

|  400 

31 

i 

1  450 

8 

i 

Table  A7 

Duration  Information  for  Newport,  OR 


Number 

Number 

H  >  HI 

H  >  HI 

of 

of 

i  cm 

% 

Events 

Events/yr 

Table  A8 


Duration  Information  for  Half-Moon  Bay,  CA 


HI 

cm 

H  >  HI 

cm 

H  >  HI 

% 

Number 

of 

Events 

500 

768 

1.3 

105 

550 

373 

0.6 

50 

600 

168 

0.3 

23 

650 

17 

0.03 

3 

Number 

of 

^min 

tmax 

t 

°t 

Events/yr 

hrs 

hrs 

hrs 

hrs 

5.2 

3 

123 

21 .7 

21  .2 

2.5 

3 

108 

22.0 

21.7 

1.2 

3 

78 

21.5 

17.8 

0.2 

12 

21 

17.0 

4.6 

A24 


Table  A9 


Duration  Parameters  with  a  Peak  Wave  Height  Threshold, 


Nagshead.  NC 


H2  %  Records 
cm  H  >  HI 


Number  Number 

of  of 

Events  Events/yr 


Table  A10 


Table  All 

Distribution  Parameters  for  Durations  and  Peak  Wave  Heights  at 


Atlantic  City,  NJ 


Parameter 


- Duration - 

Sample  Type  I  Weibull 


HI  =  200  cm 

(2.5% 

occurrence  level) 

c/a 

- 

7.75 

1  .2 

8 

- 

8.93 

14. C 

X 

12.8 

12.9 

13.  1 

0 

10.9 

11.4 

10.7 

r 

- 

0.97 

0.9 

2* 

Eres 

- 

1 .40 

1.4 

std.err. 

- 

0.066 

O.C 

HI  =  250  cm 

(0.7% 

occurrence ) 

c/a 

- 

6.36 

1.4 

B 

- 

6.23 

11.  C 

X 

9.9 

10.0 

10. C 

0 

7.3 

8.0 

7.C 

r 

- 

0.97 

0.9 

2 • 

Eres 

- 

0.512 

0.5 

3td.err . 

0.068 

O.C 

{ Continued ) 

•Sum  of  the  square  residuals 


- Peak  Wave  Height — 

Sample  Type  I  Weibull 


2/2 


”o-A184  055  DURATION  OF  EXTREME  NAVE  CONDITIONSCU)  COASTAL 

ENGINEERING  RESEARCH  CENTER  VICKSBURG  MS  0  P  SMITH 
JUN  87  CERC-MP-87-12 


UNCLASSIFIED 


F/O  8/3 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-  1963-A 


Table  All  (Concluded) 


- Duration -  - Peak  Wave  Height — 

Parameter  Sample  Type  I  Weibull  Sample  Type  I  Wei bull 

HI  s  300  cm  <0.1%  occurrence) 


e/a 

- 

4.65 

1.45 

- 

318.7 

12.6 

B 

- 

6.39 

9.20 

- 

25.2 

345.5 

X 

8.1 

8.3 

8.6 

332.2 

333.2 

331.6 

a 

6.8 

8.2 

7.1 

27.7 

32.3 

32.0 

r 

- 

0.92 

0.92 

- 

0.98 

0.97 

2* 

Eres 

- 

0.366 

0.366 

- 

0.067 

0.139 

std.err . 

- 

0.116 

0.116 

- 

0.050 

0.072 

HI  =  350  cm 

(0.03*  occurrence) 

e/a 

- 

3.43 

1.36 

- 

354.5 

15.4 

S 

- 

4.57 

6.78 

- 

21.8 

376.6 

X 

5.7 

6.1 

6.2 

365.2 

367.1 

364.0 

a 

4.1 

5.9 

4.6 

18.4 

28.0 

29.0 

r 

- 

0.89 

0.89 

- 

0.86 

0.77 

2» 

Eres 

- 

0.128 

0.122 

- 

0.159 

0.243 

std.err. 

0.135 

0.132 

_ 

0  .151 

0.186 

*  Sub  of  the  square  residuals. 


A28 


Table  A12 


Distribution  Parameters  for  Durations  and  Peak  Wave  Heights  at 


Nagshead.  NC 


(  * 

1 

> 

Parameter 

Sample 

-Duration - 

Type  I  Weibull 

- Peak  Wave  Height — 

Sample  Type  I  Weibull 

HI  =  150  cm 

(11. 9%  occurrence) 

e/o 

- 

6.19 

1.06 

- 

181.5 

4.80 

't 

i 

6 

- 

36.0 

25.0 

- 

44.5 

226.8 

X 

26.9 

27.0 

24.5 

207.0 

207.2 

207.7 

a 

39.1 

46.2 

23.2 

55.1 

57.1 

49.4 

V 

r 

- 

0.88 

0.98 

- 

0.99 

0.95 

2 

Zres 

- 

14.5 

2. 18 

- 

1.59 

6.43 

$ 

std.err . 

- 

0.14 

0.05 

- 

0.04 

0.09 

• 

HI  =  200  cm 

(5.4%  occurrence) 

e/o 

- 

6.20 

0.99 

- 

218.5 

5.55 

6 

- 

26.0 

19.6 

- 

43.2 

264.8 

'< 

■*< 

X 

21.1 

21.2 

19.7 

243.2 

243.4 

244.6 

0 

29.4 

33.4 

20.0 

51.8 

55.4 

51.0 

• 

r 

- 

0.90 

0.97 

- 

0.95 

0.89 

r, 

*> 

•  ■ 

2 

Zres 

- 

7.35 

2.53 

- 

3.50 

7.81 

std.err. 

- 

0.13 

0.07 

- 

0.09 

0.13 

1 

£ 

i 

y 

y 

1 

v 

(Continued) 
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Table  A12  (Continued) 


- Duration -  - Peak  Wave  Height — 

Parameter  Sample  Type  I  Weibull  Sample  Type  I  Weibull 

HI  =  250  cm  ( 1 .9%  occurrence) 


e/a 

- 

6.21 

1.02 

- 

272.3 

6.35 

6 

- 

22.3 

17.9 

- 

45.4 

321.4 

X 

18.9 

19.1 

17.7 

298.0 

298.5 

299.1 

a 

25.0 

28.6 

17.3 

53.9 

58.2 

55.0 

r 

- 

0.90 

0.96 

- 

0.96 

0.90 

2 

Ires 

- 

2.79 

1.04 

- 

1.05 

2.73 

std.err. 

- 

0.13 

0.08 

- 

0.08 

0.12 

HI  =  300  cm 

(0.6%  occurrence) 

e/a 

- 

6.30 

1.15 

- 

326.3 

7.03 

6 

- 

15.8 

15.5 

- 

48.0 

378.0 

X 

15.1 

15.4 

14.7 

353.0 

354.0 

353.7 

a 

17.2 

20.2 

12.8 

55.7 

61.6 

59.2 

r 

- 

0.92 

0.98 

- 

0.97 

0.92 

2 

Ires 

- 

0.91 

0.30 

- 

0.35 

0.95 

std.err. 

- 

0.11 

0.06 

• 

0.07 

0.11 

(Continued) 
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Table  A12  (Continued) 


- Duration 

- Peak  Wave  Height — - 

Parameter 

Sample  Type  I  i 

Wei bull 

Sample 

Type  I  Weibull 

HI  =  350  cm 

(0.2*  occurrence) 

e/a 

3.9 

1.16 

- 

375.0 

7.9 

0 

15.2 

12.6 

- 

48.5 

426.4 

X 

12.2  12.7 

12.0 

401.2 

402.9 

401.3 

0 

14.9  19.6 

10.4 

53.4 

62.2 

60.3 

r 

0.87 

0.96 

- 

0.96 

0.90 

2 

ires 

0.68 

0.22 

- 

0.21 

0.52 

atd.err. 

0.14 

0.08 

- 

0.08 

0.12 

HI  =  400  cm 

(0.1 0J  occurrence) 

c/a 

6.8 

0.97 

- 

428.9 

8.4 

B 

12.1 

14.3 

- 

53.5 

481.7 

X 

12.9  13.8 

14.5 

456.0 

459.8 

454.5 

a 

12.0  15.5 

15.0 

53.4 

68.6 

64.8 

r 

0.96 

0.97 

- 

0.97 

0.95 

2 

ires 

0.08 

0.06 

- 

0.05 

0.10 

std.err. 

0.08 

0.08 

0.07 

0. 10 
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Table  A12  (Concluded) 


- Duration -  - Peak  Wave  Height — 

Parameter  Sample  Type  I  Weibull  Sample  Type  I  Wei bull 

HI  =  450  cm  (0.03*  occurrence) 


e/a 

- 

4.2 

1.51 

- 

460.4 

8.6 

8 

- 

4.5 

7.52 

- 

52.9 

512.2 

X 

6.4 

6.8 

6.8 

486.0 

490.9 

484.1 

a 

4.1 

5.8 

4.6 

46.3 

67.9 

66.9 

r 

- 

0.93 

0.94 

- 

0.89 

0.84 

2 

ires 

- 

0.07 

0.06 

- 

0.10 

0.16 

std.err. 

_ 

0.11 

0.10 

_ 

0.13 

0.16 

Table  A13 


Distribution  Parameters  for  Durations  and  Peak  Wave  Heights  at 

Daytona  Beach,  FL 


Duration -  - Peak  Wave  Height — 


Parameter 

HI  3  200  cm 

Sample  Type  I 

(7.2%  occurrence) 

Weibull 

Sample 

Type  I  Weibull 

e/a 

- 

12.4 

0.96 

- 

221.3 

5.63 

0 

- 

30.0 

27.9 

- 

42.4 

266.9 

X 

29.5 

29.7 

28.4 

295.3 

245.7 

246.7 

a 

35.6 

38.5 

29.7 

51.4 

54.4 

50.7 

r 

- 

0.95 

0.98 

- 

0.97 

0.90 

ires2 

- 

3.56 

1.11 

- 

2.43 

6.52 

std.err . 

HI  =  250  cm 

0.09 

(2.3%  occurrence) 

0.05 

0.08 

0.12 

e/a 

- 

10.2 

1.01 

- 

217.0 

6.48 

0 

- 

21.6 

21.9 

- 

42.8 

318.1 

X 

22.5 

22.7 

21.8 

295.3 

295.7 

296.3 

o 

25.5 

27.6 

21.7 

51.4 

54.9 

53.5 

r 

- 

0.95 

0.98 

- 

0.97 

0.90 

2 

Ires 

- 

1.51 

0.68 

- 

1.02 

2.79 

std.err. 

_ 

0.09 

0.06 

- 

0.07 

0. 12 

(Continued) 


A33 


Table  A13  (Concluded) 


-Duration - 

- Peak  Wave  Height —  1 

Parameter 

Sample 

Type  I 

Weibull 

Sample 

Type  I  Weibull 

HI  s  300  cm 

(0.8%  occurrence) 

e/a 

- 

10.8 

1.21 

- 

329.7  8.43 

0 

- 

15.5 

20.6 

- 

39.4  372.6 

X 

19.5 

19.8 

19.3 

351.6 

352.4  351.7 

a 

18.0 

19.9 

16.1 

46.4 

50.5  49.7 

r 

- 

0.98 

0.99 

- 

0.99  0.93 

2 

Ires 

- 

0.29 

0.11 

- 

0.17  0.78 

std.err. 

- 

0.06 

0.04 

- 

0.05  0.10 

X 

u 

o 

i 

(0.2%  occurrence) 

e/a 

- 

7.39 

1.40 

- 

375.7  10.4 

0 

- 

10.9 

14.6 

- 

35.4  413.9 

X 

13.3 

13.7 

13.3 

394.8 

396.2  394.4 

o 

11.5 

14.0 

9.61 

39.4 

45.4  45.7 

r 

0.95 

0.98 

- 

0.98  0.93 

ires^ 

- 

0.27 

0.10 

- 

0.13  0.36 

std.err. 

1 

1 

i 

• 

0.09 

0.06 

0.06  0.11 

a 
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Table  Al4 


Distribution  Parameters  for  Durations  and  Peak  Wave  Heights  at 

Newport,  OR 


Duration -  - Peak  Wave  Height — 


Parameter  Sample  Type  I  Welbull  Sample  Type  I  Weibull 

HI  =  400  cm  (17. 9%  occurrence) 


e/a 

- 

15.4 

1.03 

- 

457.9 

7.45 

6 

- 

39.7 

36.2 

- 

62.7 

527.3 

X 

38.2 

38.4 

35.7 

493.9 

494.1 

494.8 

a 

47.1 

51.0 

34.6 

77.7 

80.4 

78.5 

r 

- 

0.94 

0.99 

- 

0.99 

0.95 

2 

ires 

- 

8.07 

0.87 

- 

1.52 

6.81 

std.err . 

- 

0.10 

0.03 

- 

0.04 

0.09 

HI  =  450  cm 

(11.1%  occurrence) 

e/a 

- 

13.6 

1.08 

- 

505.5 

7.51 

B 

- 

28.6 

29.4 

- 

53.5 

565.4 

X 

30.1 

30.2 

28.6 

536.2 

536.4 

536.7 

a 

34.3 

36.7 

26.5 

65.9 

63.6 

67.7 

r 

- 

0.96 

0.99 

- 

0.99 

0.96 

2 

Eres 

- 

4.75 

0.92 

- 

1.55 

3.96 

std.err. 

0.09 

0.04 

- 

0.05 

0.08 

(Continued) 
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Table  A14  (.Continued) 


Parameter 

Sample 

-Duration 
Type  I 

Weibull 

- Peak  Wave  Height — 

Sample  Type  I  Weibull 

HI  =  500  cm 

(3.7)1  occurrence) 

e/a 

- 

11.8 

1.13 

- 

547.4 

13.0 

8 

- 

22.3 

24.7 

- 

41.6 

594.6 

s- 

X 

24.5 

24.6 

23.6 

571.3 

571.5 

571.5 

a 

26.5 

28.6 

21.0 

51.4 

53.4 

53.5 

r 

- 

0.96 

0.99 

- 

0.98 

0.98 

t  p 

2 

Ires 

- 

3.01 

0.72 

- 

1.24 

1.68 

Y 

std.err. 

- 

0.08 

0.04 

- 

0.05 

0.06 

f 

.  < 

HI  =  550  cm 

(3.7%  occurrence) 

e/a 

- 

10.8 

1.21 

- 

585.0 

18.9 

0 

- 

14.9 

20.2 

- 

29.8 

619.3 

X 

19.3 

19.4 

18.9 

602.0 

602.2 

602.0 

t 

o 

18.2 

19.2 

15.7 

37.1 

38.2 

39.5 

r 

- 

0.98 

0.99 

- 

0.99 

0.97 

**, 

Ires2 

- 

1.32 

0.61 

- 

0.40 

1.47 

■Y 

■  t 

Y\ 

V, 

c  V 

' 

std.err. 

0.06 

0.04 

0.03 

0.07  < 

v 

i1'. 
1* 
t*  ; 

(Continued) 
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Table  Al4  (Continued) 


- Duration -  - Peak  Wave  Height — 

Parameter  Sample  Type  I  Wei bull  Sample  Type  I  Weibull 

HI  =  600  cm  (1.8)1  occurrence) 


e/a 

- 

9.1 

1.13 

- 

616.6 

24.4 

6 

- 

13.1 

17.2 

- 

22.6 

643.6 

X 

16.5 

16.7 

16.5 

629.5 

629.7 

629.5 

a 

15.8 

16.8 

14.6 

27.6 

29.0 

32.1 

r 

- 

0.97 

0.98 

- 

0.99 

0.98 

2 

Ires 

- 

1.09 

0.75 

- 

0.41 

2.53 

std.err . 

- 

0.07 

0.06 

- 

0.05 

0.11 

HI  =  650 

cm  (0.3*  occurrence) 

e/a 

- 

5.1 

1.12 

- 

663.7 

32.2 

B 

- 

10.1 

11.4 

- 

18.8 

685.1 

X 

10.6 

11.0 

10.8 

673.9 

674.5 

673.5 

a 

10.9 

13.0 

9.2 

21.4 

24.1 

26.3 

r 

- 

0.92 

0.95 

- 

0.98 

0.93 

2 

Ires 

- 

0.52 

0.34 

- 

0.13 

0.50 

std.err. 

0.11 

0.09 

- 

0.06 

0.11 

(Continued) 
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Table  A14  (Concluded) 


- Duration -  - Peak  Wave  Height — 

Parameter  Sample  Type  I  Weibull  Sample  Type  I  Wei bull 

HI  =  700  (0.0451  occurrence) 


e/o 

- 

5.4 

1.15 

- 

711.5 

71.9 

8 

- 

9.4 

11.5 

- 

10.4 

721.1 

X 

9.9 

10.8 

10.9 

716.4 

717.5 

715.5 

o 

8.3 

12.1 

9.5 

9.2 

13.3 

12.6 

r 

- 

0.93 

0.96 

- 

0.97 

0.98 

2 

tree 

- 

0.06 

0.03 

- 

0.03 

0.02 

std.err. 

0.11 

0.08 

0.07 

0.06 

Table  A15 

Distribution  Parameters  for  Durations  and  Peak  Wave  Heights  at 

Half-Moon  Bav,  CA 


Duration -  - Peak  Wave  Height — 


Parameter  Sample  Type  I  Weibull  Sample  Type  I  Weibull 

HI  =  500  cm  (1.3/t  occurrence) 


e/a 

- 

11.7 

1.11 

- 

533.1 

14.2 

6 

- 

18.0 

22.5 

- 

36.3 

574.0 

X 

21.7 

22.0 

21.7 

553.4 

554.0 

553.3 

0 

21 .2 

23.1 

19.6 

43.7 

46.6 

47.6 

r 

- 

0.98 

0.99 

- 

0.99 

0.97 

2 

ires 

- 

0.42 

0.19 

- 

0. 12 

0.58 

std.err . 

- 

0.06 

0.04 

- 

0.03 

0.08 

HI  =550  cm 

(0.6/1  occurrence) 

e/a 

- 

11.4 

1.14 

- 

575.9 

19.4 

3 

- 

19.4 

22.9 

- 

28.3 

607.4 

X 

22.0 

22.6 

21.8 

591.3 

592.1 

590.8 

0 

21.7 

24.9 

19.2 

32.4 

36.1 

37.7 

r 

- 

0.96 

0.99 

- 

0.99 

0.95 

2 

Eres 

- 

0.32 

0.09 

- 

0.06 

0.39 

std.err. 

0.08 

0.04 

- 

0.04 

0.09 

Table  A15  (Concluded) 


Parameter  Sample  Tyc 

HI  =  600  cm  (0.3J  occurrence) 


- Duration - 

Sample  Type  I  Wei bull 


- Peak  Wave  Height — 

Sample  Type  I  Welbull 


e/a 

- 

12.8 

1.17 

- 

607.1 

20.0 

8 

- 

16.6 

23.7 

- 

26.9 

637.4 

X 

21.5 

22.3 

22.4 

621.3 

622.7 

620.5 

a 

17.8 

21.3 

19.3 

28.0 

34.5 

38.4 

r 

- 

0.99 

0.99 

- 

0.93 

0.86 

2 

£res 

- 

0.05 

0.02 

- 

0.22 

0.46 

std.err. 

- 

0.05 

0.03 

- 

0.10 

0.15 

r. 

n 


f. 


i 
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Table  A1 6 


Table  Al8 


i 


Results  of  Regression  of  Duration  Against  Conditions  at 
the  Peak  of  the  Event  for  Daytona  Beach.  FL 


2 

2 

2  2 

2 

HI 

wBSBSm, 

rH 

rH 

rT 

rT 

rH,H  tT,T 

rH  L 

200 

7.1 

0.76 

0.73 

0.50 

0.48 

0.79 

0.65 

250 

2.2 

0.81 

0.79 

0.46 

0.46 

0.82 

0.72 

300 

0.8 

0.51 

0.50 

0.46 

0.47 

0.59 

0.52 

350 

0.2 

0.44 

0.46 

0.52 

0.55 

0.66 

0.59 

i 

I 


Table  A19 

Results  of  Regression  of  Duration  Against  Conditions  at 
the  Peak  of  the  Event  for  Newport .  OR 


Table  A20 


Results  of  Regression  of  Duration  Against  Conditions  at 


the  Peak  of  the  Event  for  Half-Moon  Bay.  CA 


2 

1  rH 

rH 

r1 

’  2  2 

2 

r  rH,H  ,T,T 

rH  L 

APPENDIX  B 


PERTINENT  DATA  FROM  THE  WAVE  INFORMATION  STUDIES  PROGRAM 


Figures  Page 


B1  Phase  III  wave  rose,  Atlantic  City,  NJ .  B2 

B2  Phase  III  wave  rose,  Nagshead,  NC .  B3 

B3  Phase  III  wave  rose,  Daytona  Beach,  FL .  B4 


Tables 

B1  Phase  III  Wave  Data,  Atlantic  City,  NJ: .  B2 

B2  Phase  III  Wave  Data,  Nagshead,  NC .  B3 

B3  Phase  III  Wave  Data,  Daytona  Beach,  FL .  B4 

B4  Phase  III  Wave  Data,  Newsport,  OR .  B5 

B5  Phase  III  Wave  Data,  Half-Moon  Bay,  CA .  B5 

B6  Phase  III  Duration  Data,  Atlantic  Coast .  B6 


I 

I 

I 


I 


HKXSHTIHCTRCS) 


Table  B1 

Phase  III  Wave  Data.  Atlantic  Cltv.  MJ 
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Figure  B1 .  Phase  III  wave  rose,  Atlantic  City,  NJ 
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Table  52 


Phase  III  Wave  Data.  Nagshead,  NC 
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Figure  B2.  Phase  III  wave  rose,  Nagshead,  NC 
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Phase  III  Wave  Data,  Daytona  Beach.  FL 
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Figure  B3.  Phase  III  wave  rose,  Daytona  Beach,  FL 
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Table  B4 


Phase  III  Wave  Data.  Newport,  OR 


PERCENT  OCCURRENCEc X100 ) 


HEIGHT) METERS) 
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Table  B5 

Phase  III  Wave  Data,  Half-Moon  Bay,  CA 
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FORTRAN  Listing  -  Program  "STRMDIST" 


Mi 


10**N,  J 

20*: IUENT:R0CDOPS,OPSMITH 
30*:OPT1ON:FORTRAN 
401: USE: . STL  I T 
50*: FORTY 

60C  HONEYWELL  VERSION  1/2/86 

70C  *#  PROGRAM  "STRMDIST"  READS  A  WIS  PHASE  III  DATA  FILE  ** 

80C  #*  AND  IDENTIFIES  STORMS  WHERE  CONSECUTIVE  RECORDS  ** 

90C  **  HAVE  WAVE  HEIGHTS  EXCEEDING  A  SPECIFIED  THRESHOLD.  ♦* 

100C  **  THE  NUMBER,  PEAK  CONDITIONS  AND  DURATIONS  OF  THESE  ♦♦ 

110C  **  STORMS  ARE  THEN  TABULATED.  STORMS  ONLY  6  HOURS  ** 

120C  *#  APART  ARE  CONSIDERED  AS  A  SINGLE  EVENT.  THE  ** 

130C  #t  PROGRAM  ALSO  FITS  AN  EXTREMAL  TYPE  I  AND  A  WE  I  BULL  ** 

140C  **  DISTRIBUTION  TO  THE  PEAK  WAVE  HEIGHTS  AND  THE  DUR-  ** 

150C  **  ATI ONS  AND  REPORTS  THE  PARAMETERS  OF  EACH.  t* 

1 60C 

170  DIMENSION  DUR ( 999 ) , HPE AK ( 999 ) , TPE AK ( 999 ) ,DPEAK(9°9> ,DTPEAK<999> 

180  INTEGER  DATM , HSEA , TSEA , DSEA,HSWL , TSWL , DSWL , STNO , DAT IM , H , T , D 

190  INTEGER  STMNO , ON, HPK , TPK , DPK , HPEAK , TPEAK , DPEAK, DUR , HI , H2 , RECNO 

200  INTEGER  DTPEAK , DAT  I  ME , DTPK , FLAG , NOREC , NOYRS , TM IN , TMAX 

210  INTEGER  YR , YRP , MO , MOP , DY , DYP , TM , TMP , YRD IFF , MODI FF , DYD I FF , 

220  VTMDIFF ,HPKMIN,HPKMAX 
230  CHARACTER*64  FNAME 

240  CHARACTER*8  VARIABLE 

250  CHARACTER*8  VARIABLE 

260C 

270C  #*  READ  WIS  DATA  FILE  AND  WRITE  FILE  OF  STORMS  ** 

280C  ♦*  EXCEEDING  1ST  WAVE  HEIGHT  THRESHOLD,  HI  (CM)  ** 

290C 

300  FNAME=  “ NAGSHEAD ,  NORTH  CAROLINA" 

310  NOREC=58440 

320  NOYRS=20 

330  H 1=300 

340C  **  K  »  THE  STORM  NO.  ASSIGNED  TO  CONSECUTIVE  RECORDS  ** 

350  K=0 

360C  **  3  »  THE  NO.  RECORDS  WHERE  H  >  HI  ** 

370  J=0 

380  YRP=999 


.4 


% 


a 


MQP=999 
DYP=999 
TMP  =999 

CALL  ATTACH! 03 , "/A  3  083; ",  1 ,0, 1  ST AT) 

ISTAT=FLD<6,6, ISTAT) 

IF  ( ISTAT.NE.0)  GO  TO  900 
CALL  FMEDI A (07,6) 

CALL  FMEDI  A  ( 08 , 6 ) 

CALL  FMED I A ( 09 , fc ) 

READ  ll, 10)  YR,MO,iM  ,  TM,HSEA,TSEA,DSEA,HSWL, TSWL, DSWL 
FORMAT ( 2  X ,4I2,6lb/ 


i 


a 

.4 


580 C  •*  COMPUTE  COMPOSITE  SEA  AND  SWELL  WAVE  HEIGHT,  1ST  RECORD  *# 

510C 

528  H»INT (SORT (FLOAT (HSEA)**2+FL0AT(HSWL>**2) ) 

538  IF(H.LT.Hl)  60  TO  38 

548  J-l 

558  K-l 

568C 

S78C  ••  SET  COMPOSITE  PERIOD  AND  DIRECTION  TO  THAT  OF  SEA  OR  ** 

S88C  ••  SWELL,  WHICHEVER  HAS  A  HIGHER  INCREMENTAL  WAVE  HEI6HT  ## 

59BC 

688  IF(HSEA.GT.HSNL)  GO  TO  15 

618  T-TSNL 

628  D-DSWL 

638  GO  TO  20 

640  15  T-TSEA 

650  D-DSEA 

660  20  WRITE (7,25)  K , YR , MO , DY , TM, H, T , D 

670  25  FORMAT  <2X, 14, IX, 412, 316) 

680  YRP-YR 

698  MOP-MO 

708  DYP-DY 

710  TMP-TM 

728  30  READ(1, 10, END-200)  YR , MO , OY , TM , HSEA , TSEA , DSE A , HSWL , TSWL , DSWL 

730C 

740C  **  COMPUTE  COMPOSITE  SEA  AND  SWELL  WAVE  HEIGHT  *• 

750C 

760  H-INT (SORT (FLOAT (HSEA) **2+FL OAT (HSWL) **2) ) 

770  32  IF(H.LT.Hl)  GO  TO  190 

780  J-J  +  l 

790C 

B00C  **  SET  COMPOSITE  PERIOD  AND  DIRECTION  TO  THAT  OF  SEA  OR  ** 

810C  *#  SWELL,  WHICHEVER  HAS  A  HIGHER  INCREMENTAL  WAVE  HEIGHT  ## 

B20C 

830  IF (HSEA. GT. HSWL)  GO  TO  35 

840  T -TSWL 

850  D-DSWL 

860  GO  TO  40 

870  35  T-TSEA 

880  D-DSEA 

890  40  YRDIFF-YR-YRP 

900  MODIFF-MOP-MO 

910  DYDIFF-DYP-DY 

920  TMDIFF-TMP-TM 

930C 

940C  ***  CHECK  FOR  CONSECUTIVE  RECORDS  (SAME  STORM)  *#* 

950C 

960C  **  CONSECUTIVE  RECORDS,  SAME  DAY  ** 

970  IF(YRDIFF.EQ.0. AND . MODI FF . EQ . 0 . AND . DYD I FF . EQ . 0 . AND. 

980  6TMDIFF. EQ. -3)  GO  TO  45 

990C  **  CONSECUTIVE  RECORDS,  DAY  END  *♦ 

1000  I F ( YRDIFF . EQ. 0. AND. MOD  IFF . EQ. 0. AND. D YD  I FF . EQ. - 1 . AND. 

1010  6TMD IFF . EQ . 2 1 )  GO  TO  45 

1020C  **  CONSECUTIVE  RECORDS,  MONTH  END  *« 


1039  IF(YRDIFF. EQ.0. AND. HOOIFF. EQ.-l. AND. OYDIFF.EQ, 27. AND 

1040  UMDIFF.EO.2l)  80  TO  43 

1030  IF ( YRDIFF. EQ.O. AND. HODIFF. EQ.-l. AND. OYDIFF.EQ. 28. AND 

1860  6THDIFF. EQ. 21 )  60  TO  43 

1070  IF (YRDIFF. EQ.O. ANO. HODIFF. EQ.-l. AND. DYDIFF.EQ. 29. AND. 

1080  6TflD1FF.EQ.21)  60  TO  43 

1090  IF (0ROIFF.EQ.O. ANO. HOOIFF. EQ. -1. AND. OYDIFF.EQ. 38. ANO. 

1100  6TN01FF.EQ.21)  60  TO  43 

1 1 10C  **  CONSECUTIVE  RECORDS,  YEAR  END  #• 

1120  IF ( YRDIFF. EQ.l. AND. flOOIFF.EQ.il. AND. OYDIFF.EQ. 30. AND. 

1130  6TN01FF.EQ.21)  30  TO  43 

1140C 

1130C  ***  CHECK  FOR  RECORDS  6  HRS  APART  AND  ADJUST  RECORD 

116BC  ***  BETWEEN  SUCH  THAT  THE  PR06RAN  SEES  ONE  CQNT- 

1170C  INUOUS  STORN  (I6N0RIN8  THE  ONE  RECORD  BELOW 

1180C  ***  THE  THRESHOLD) 

1 190C 

1208C  *•  RECORDS  6  HRS  APART,  SANE  DAY  #• 

1218  IF ( YRDIFF. EQ.0. AND. HODIFF.EQ.0. AND. OYDIFF.EQ.0. AND. 

1220  6TNDIFF.EQ.~6)  60  TO  47 

1230C  *•  RECORDS  6  HRS  APART,  DAY  END  *t 

1240  IF (YRDIFF. EQ.0. AND. HODIFF.EQ.0. AND. DYDIFF.EQ. -1. AND. 

1230  6TNDIFF.EQ. 10)  80  TO  47 

1260C  •»  RECORDS  6  HRS  APART,  HONTH  END  •• 

1278  IF (YRDIFF. EQ.0. AND. H001FF.EQ.~1. AND. OYDIFF.EQ. 27. AND. 

1200  6TH0IFF . EQ. 18)  60  TO  47 

1290  IF (YRDIFF. EQ.0. AND. NOOIFF.EQ.-l. AND. OYDIFF.EQ. 20. AND. 

1300  6TNDIFF. EQ. 18)  80  TO  47 

1310  IF ( YRO I FF. EQ. 0. AND. HOD IFF. EQ.-l. AND. OYDIFF.EQ. 29. ANO. 

1320  6THDIFF.EQ. 18)  60  TO  47 

1330  IF (YRDIFF. EQ.0. ANO. NOOIFF. EQ.-l. AND. OYDIFF.EQ. 30. AND. 

1340  6TNDIFF. EQ. IS)  60  TO  47 

I330C  *•  RECORDS  6  HRS  APART,  YEAR  END  •• 

1360  IF (YRDIFF. EQ.l. ANO. NODIFF.EQ.il. ANO. OYDIFF.EQ. 30. ANO. 

1370  6TNDIFF.EQ. IS)  60  TO  47 

1380  K-KM 

1390  60  TO  43 

1400  47  BACKSPACE  1 

1410  BACKSPACE  1 

1420  READ (1,10)  YR ,RO,DY ,TH,HSEA,TSEA,DSEA,HSNL,TSWL,DSHL 

1430  H»H1 

1440  60  TO  32 

1431  43  Nft!TE(7,23)  K, YR,NO,DY ,TN,H,T ,0 

1460  YRP-YR 

1470  HOP«HO 

1401  DYP«OY 

1490  THP«TH 

1301  190  80  TO  30 
1310C 

1320C  ••  FILE  CODE  7  INCLUDES  RECORDS  WHERE  H  IS  6REATER  THAN 

1330C  ••  THE  FIRST  WAVE  HEI6HT  THRESHOLD  HI.  CONSECUTIVE 

1340C  **  RECORDS  SHARE  A  COHHON  ‘STORH  NUMBER*,  K. 


|  1S60C 

I  157BC 

*«  READ  FILE  OF  STORMS  ,  COMPUTE  DURATIONS  AND 

*# 

■  158BC.  t*  IDENTIFY  PEAK  CONDITIONS 

•  • 

I  1398C 

[  1610 

208  RECNO-J 

1  1610 

REMIND  7 

j  1^20 

DN»0 

|  1630 

HPK*0 

1  1640 

J«1 

I  1650 

95  READ (7, 97, END" 400)  STNO,DATH,H,T,D 

I  1660 

97  FORMAT <2X, 14, IX, 18,316) 

B  1670 

IF(STNO.EQ.J)  60  TO  3S0 

1  1600 

DUR ( J ) *DN 

R  1690 

DTPEAK ( J ) «DTPK 

I  1700 

HPEAK ( J ) *HPK 

I  1710 

TPEAK ( J ) "TPK 

1  1720 

DPEAK ( J) "DPK 

1730 

J-J  +  l 

1740 

DN«0 

1750 

HPK-0 

1760 

350  DN-DN+3 

1770 

IF(H.LE.HPK)  GO  TO  95 

I  1780 

DTPK-DATM 

1790 

HPK-H 

*  1000 

TPK«T 

\  1810 

DPK*D 

j  1820 

GO  TO  95 

1830 

400  DUR ( J) *DN 

1840 

DTPEAK  <  J ) »DTPK 

1850 

HPEAK ( J ) “HPK 

1860 

TPEAK ( J ) "TPK 

1870 

DPEAK ( J ) *DPK 

1880 

PQISSON-FLOATtJ) /FLOAT (NOYRS) 

1890 

PERCENT«RECNO*10B. /NOREC 

1900 

TMIN-DUR  < 1 > 

i  1910 

HPKHIN-HPEAK(l) 

1920 

HPKHAX*H1 

1930 

TMAX-3 

1940 

TSUH»0.0 

1950 

HPKSUM-0.0 

1960 

DO  700  I«1,J 

1970 

TSUH»TSUH+FLOAT l  DUR  !  I ) ) 

1980 

HPKSUM-HPKSUM+FLOAT (HPEAK ( I ) > 

1990 

IF (HPEAK (I).LT.HPKMIN)  HPKMIN»HPEAK ( I ) 

2000 

IF ( HPEAK (I).BT.HPKMAX)  HPKHAX=HPEAK l I ) 

2010 

IF (DUR ( I ) . LT. THIN)  TMIN«OUR(I) 

n 

2020 

IF(DUR(I) .GT.TMAX)  TMAX»DUR(I) 

2030 

700  CONTINUE 

1 

2040 

HPKHEAN»HPKSUH/F10AT (J) 

2050 

TMEAN»TSUM/FLOAT(J) 

2060 

TDIFFSUM«0.0 

2070 

HDIFSUM-0.0 

2080 

DO  710  I • 1 , J 

C5 

: 

j 

rvrvr 


1*,  6  r 


.;^w 


2891  HD1FSQ- (FLOAT (HPEAK  < I ) ) -HPKMEAN) **2 

2188  HDIFSUM-HDIFSUM+HDIFSQ 

2118  TDIFSfl* (FLOAT (OUR ( I ) ) -TMEAN) «#2 

2128  TDIFFSUM"TDIFFSUH+TDIFSQ 

2138  718  CONTINUE 

2148  STDEVT«SBRT (TDIFFSUH/FLOAT  <  J - 1 ) ) 

2138  STDEVH"SBRT(HD1FSUM/FL0AT<J-1) ) 

2168C 

217BC  •#  PRINT  TABLE  OF  STORN  PARAMETERS  •• 

21  see 

2198  NRITE (6,448)  FNAME 

2288  448  FORMAT (1H1 ,///, 25X , "ANALYSIS  OF  STORM  DURATION", 

2218  fc//,8X ,  "DATA  FILE:  “,A64) 

2228  WRITE (6,438)  HI 

2238  438  FORMAT (//, IX , “STORM  NO. " , 2X , "DATE/T IME  OF  PEAK  ",2X, 

2248  INDURATION  H>" , 13, 2X , “PEAK  H“,2X,"PEAK  T",2X,"PEAK  DIR",/) 

2258  DO  388  L- 1 , J 

2268  WRITE (6,478)  L,DTPEAK(L) ,DUR(L) , HPEAK (L) ,TPEAK (L) , DPEAK (L) 

2278  WRITE (9 , 470)  L , DTPEAK (L ), DUR <L> , HPEAK (L) , TPEAK (L) , DPEAK (L ) 

2288  478  FORMAT <4X , 1 3 , 10X , 18 , 0X , 14 , 12X , 1 3 , 6X , 1 2 , 6X , 1 3) 

2298  388  CONTINUE 

2388  WRITE (6,510)  POISSON 

2318  518  FORMAT ( /  , 4X  ,F3. 2 , "  STORMS  PER  YEAR") 

2328  NRITE (6,28)  HI ,RECNO, PERCENT, NOREC 

2338  28  FORMAT (/ ,4X , "NO.  RECORDS  WHERE  H  >  ",I3,“  «  ",I5, 

2340  i"  ( " ,F4. 1 , "X  OF  *,I5,"  RECORDS)") 

2358  WRITE (6,719)  HPKH IN , HPKMAX , HPKMEAN , STDEVH 

2368  719  FORMAT ( / , 4X , "MIN.  PEAK  H  -  “ , 1 3 , -  MAX.  «  ",I4, 

2378  I."  MEAN  ■  ",F5.1,"  STD.  DEV.  -  ",F5.1> 

2380  WRITE (6,720)  THIN, TMAX .TMEAN, STDEVT 

2398  728  FORMAT (/, 4X , "MIN.  DURATION  -  " ,  13 , "  MAX.  «  “,I3, 

2488  l"  MEAN  -  “,F5.1,“  STD.  DEV.  »  ",F5.1> 

2418  WRI TE (6 , 475) 

2420  475  FORMAT (//,4X,“THE  DATE/TIME  IS  YRMODYHR ,  DURATION  IS  IN  HOURS, 

2430  liH  (HEIGHT)  IS  IN  CM,",/,4X,"T  (PERIOD)  IS  IN  SEC  AND  DIRECTION" 
2448  6"  IS  IN  DEGREES  RELATIVE  TO  THE  SHORELINE") 

2438  VARIABLE-'PEAK  H  ' 

2468  CALL  PROBDIST  (VARIABLE,.], HPEAK, POISSON) 

2478  VARIABLE" 'OURAT I ON ' 

2488  CALL  PROBDIST  (VARIABLE, J , DUR, POISSON) 

2498  60  TO  628 

2588  988  PRINT  901 

2518  981  FORMAT ( IX , 19HATTACH  UNSUCCESSFUL) 

2520  CALL  DETACH(01 , , > 

2538  620  STOP 

2548  END 

2358C 

2368C 

237BC 

2388  SUBROUTINE  PROBDIST  ( VARIABLE ,N,HS, LAMBDA) 

2398C  SUBROUTINE  PROBDIST  ADAPTED  1/86  BY  ORSON  P.  SMITH  FROM 
2688C  PROGRAM  "WAVD1ST1".  11/85  VERSION  BY  ROBERT  B.  LUND 
2610C  DESIGN  BRANCH-COASTAL  ENGINEERING  RESEARCH  CENTER 


2620C  U.S.  ARMY  EN81NEERS  WATERWAYS  EXPERIMENT  STATION 

2630C  P.O.  BOX  631 

2640C  VICKSBURG,  MS  39181-0631 

26S0C  FOR  FURTHER  INFORMATION  CONCERNING  THE  APPLICATION 

266IC  OF  "NAVDIST1" ,  CALL _ 

267BC  ROBERT  B.  LUND  (601 ) -634-2068  FTS:206B 
2600C  ORSON  P.  SMITH  (601 ) -634-201 3  FTSt  542-2013 
2690C  DOYLE  L.  JONES  (601 ) -634-2069  FTSi542-2069 

2700C 

2710C  FORTRAN  4  HONEYWELL  DPS-8 

2720C  REFi  "RELIABILITY  OF  LONG-TERM  NAVE  CONDITIONS  PREDICTED  WITH  DATA  SETS 
2730C  OF  SHORT  DURATION"  CETN-I-5 

2740C  REFs  "HANDBOOK  OF  MATHEMATICAL  FUNCTIONS"  BY  ABRAMOWITZ  ANO  SE6UN 
27S0C  REF:  "EXTREMAL  PREDICTION  IN  NAVE  CL IMAT0L06Y"  BY  B0R6MAN  AND  RESIO 
2760C  REF.  "LONG-TERM  DISTRIBUTIONS  OF  OCEAN  WAVES." 

2770C  ISAACSON  AND  MACKENZIE 

2780C 

2790C  N  -  NUMBER  OF  STORMS 

2800C  RET  >  RETURN  PERIOD 

2810C  LAMBDA  >  POISSON  LAMBDA  PARAMETER  ( AVERA6E  NO.  STORMS  PER  YEAR) 

2820C  HS  -  THE  INDEPENDENT  VARIABLE 

2830C  DIFF  ■  THE  RESIDUAL  FOR  EACH  DATA  POINT 

2B40C  YACT  *  THE  PROBABILITY  AS  ESTIMATED  BY  THE  PLOTTING  FORMULA  M/K  +  l 

28S0C  YEST  -  THE  PROBABILITY  AS  ESTIMATED  BY  THE  DISTRIBUTION 

2860C  ALPHA  -  THE  ARRAY  OF  LOCATION  PARAMETERS  FOR  THE  DISTRIBUTIONS 

2870C  BETA  *  THE  ARRAY  OF  SCALE  PARAMETERS  FOR  THE  DISTRIBUTIONS 

28B0C  A  -  THE  SLOPE  OF  EACH  “PLOTTED  LINE" 

2890C  B  «  THE  Y-INTERCEPT  OF  EACH  “PLOTTED  LINE" 

2900C  C  »  THE  ARRAY  OF  COEFFICIENTS  FOR  THE  GAMMA  INTEGRAL  EXPANSION 

2910C  ST  «  THE  SUM  OF  THE  SQUARE  RESIDUALS 

2920C  CORR  »  THE  NON-LINEAR  CORRELATION  FOR  EACH  DISTRIBUTION 

2930C  STE  ■  THE  STANDARD  ERROR  OF  THE  ESTIMATE  OF  Y  DN  X 

2940C  MSD  »  THE  MEAN  SQUARE  DEVIATION 

29S0 

2960C  DECLARATION  OF  VARIABLES,  FUNCTIONS,  AND  CHARACTERS 

2970  DIMENSION  YACT (999,3) , YEST (999,3) , DUM1 (999) ,DUM2(999) ,HS(999) 

29B0  DIMENSION  YAVG ( 3) , CORR (3) , ALPHA ( 4) , BETA ( 4) , VAR ( 4) , DM (3 ) 

2990  DIMENSION  RET ( 5) , CHS (5 , 3) , A ( 3) , B ( 3) ,ST (3) , SB ( 3) , STE ( 3) 

3000  DIMENSION  STDEV<3) 

3010  REAL  MEAN ( 3 ) , MSD  ( 3 ) 

3020  REAL  LAMBDA 

3030  INTE6ER  HS 

3040 

3050  FI (X)»EXP(-EXP(- (X-EPSI)  /PHI)  ) 

3060  F2(X)»1.0-EXP( < - ( X/SIGMA) **C ) ) 

3070  F3(X)»£XP(-( tSI6MA2/X)**U)) 

3080 

3090  CHARACTER»20  I  FLAG  (  4 ) 

3100  CHARACTER* 1 7  DEF 

3110  CHARACTERS  F0RM<3) 

3120  CHARACTER*24  TITLE 

3130  CHARACTER* 1  L06IC 

3140  CHARACTERS  BOX  <  1 6 ) 


3150  CHARACTERS  VARIABLE 

3160C  INITIALIZATION  OF  STRINSS  AND  CONSTANTS 
3171  IFLAS < 1 ) ■ 'EXTREMAL  TYPE  I 

3180  IFLA8 (2) * ' WEI  BULL ' 

3190  1FLA6 (3) ■ ‘ LOS  EXTREMAL  ' 

3200  DEF* 'F (x ) -Pr <  X<x )  ■  ' 

3210  FORM ( 1 ) ■ ‘ EXP (-EXP ( - (x-EPSI ) /PHI ) ) 

3220  FORM (2) -‘ 1-EXP <- <x/BETA) #*ALPHA) ' 

3230  FORM (3) • ' EXP  < - (BETA/x ) t* ALPHA)  ' 

3240  TITLE-LEAST  SQUARES  RESULTS  -  ' 

3230 

3260  DATA  RET  /3. 0 , 10. 0 , 23. 0 ,50. 0 , 1 00. 0/ 

3270  EULER-. 5772156649 

3280  C2-. 7796968 

3290 
3300C 

3310C  **  SET  LOGIC  *  ‘V  FOR  PRINTOUT  OF  RESIDUAL  TABLES  «• 

3320C 

3330  LOGIC- ‘ N  ‘ 

3340C  RANK  DATA  AND  ASSIGN  A  PROB.  OF  NON-EXCEEDENCE  TO  EACH 
3350  CALL  ORDER (HS,N) 

3360  DO  25  I-l.N 

3370  DO  25  K>1,3 

3380  YACT < I, K)«FLOAT(I> /FLOATIN+l) 

3390  25  CONTINUE 

3400 

3410C  INITIALIZE  VARIABLES  FOR  LEAST  SQUARES  FIT  OF  THE  DISTRIBUTIONS 

3420  SX-0 

3430  SY-0 

3440  SXX-0 

3430  SLX-0 

3460  SLLY-0 

3470  SLXX-0 

3480  SLLQV-0 

3490  SXLLY-0 

3300  SLXLLY-0 

3510  TOOB16-0 

3520 

3330C  CALCULATE  SUMS  FOR  THE  LEAST  SQUARES  METHOD 
3540  DO  40  J-l ,N 

3550  SX-SX+HS ( J ) 

3360  SY-SY+YACT ( J  ,  1 ) 

3570  SXX»SXX*HS(J>**2 

3380  SLX-SLX*ALOS ( HS ( J ) ) 

3390  SLX  X-SLXX+ ( ALOG (HS ( J ) ) )*#2 

3600  SLLY-SLLY-ALOG ( -ALOG ( YACT  ( J  ,  1 ) ) ) 

3610  SLL0Y"SLLQY+ALOG(-ALOG(1.0-YACT (J,l) ) ) 

3620  SXLLY-SXLLY-HS(J) * ALOG ( -ALOG  (YACT  < J  ,  t ) ) ) 

3630  SLXLLY-SLXLLY- ALOB<HS<<D  ) *ALOG l -ALOG ( YACT  l  J  ,  1 ) ) ) 

3640  40  TOOBIB-TOOBIG+’ALOG(HS(J)  >*(ALOG<-  ALOG  (1.0-YACT  (J,  li  )  )  ) 

3630 

3660C  CALCULATE  SLOPE  AND  INTERCEPT  OF  EACH  "PLOTTED  LINE’ 

3670  A<1)-(N*SXLLY-SX*SLLY)/ (N#SXX-SX#*2> 


3681  AI2)»<N#TQ0BIS-SLX*SLLQY)  /  lN*SLXX-SU*»2> 

369BC  A(3)»(N*SLXLLY-SLX*SLLY> / (N*SLX X-SLX*«2) 

3711  B(1>-(SXX*SLLY-SXLLV*SX)  /  ( N * S X X - S X ■* * 2 ) 

3718  B(2)«(SLXX*SLLQY-T00BIG*SLX> / <N*SLXX-SLX*#2) 

372IC  B(3)-(SLXX*SLLY-SLXLLY*SLX)/(N«SLXX-SLX**2> 

3738C  CALCULATE  PARAMETERS  OF  EACH  DISTRIBUTION  FROM  SLOPE  AND  INTERCEPT  DATA 
3741  PH!>1.0/A{1> 

3738  EPSI*-B( 1 > / A  < 1 ) 

3768  C- A ( 2 ) 

3778  SI8MA«EXP(-B(2) /A (2) ) 

3788C  U*A ( 3 ) 

3798C  S I GMA2=EXP  <-B(3)/A(3>) 

3888 

3818C  ASSIGN  ARRAYS  ALPHA  AND  BETA  THE  PARAMETERS  OF  EACH  DISTRIBUTION 

3828C  FOR  EASY  PRINTOUT  OF  DATA 

3838  ALPHA ( 1 ) =  EPS  I 

3848  BETA ( 1 ) “PHI 

3838  ALPHA<2) *C 

3868  BETA (2) “SIGMA 

3878C  ALPHA ( 3  >  *U 

3888C  BET  A ( 3 ) -SIGMA2 

3898C  CALCULATE  PROBABILITY  AS  ESTIMATED  BY  DISTRIBUTION 

3988  DO  180  J*1,N 

3911  YEST ( J , 1 ) *F 1 ( HS  »  J  >  ) 

3928  YEST(J,2)«F2(HS(J> > 

3938C  YEST ( J , 3) *F3  <HS ( J ) ) 

3948  188  CONTINUE 

3938 

3968C  CALCULATE  AVERAGE  PROBABILITY  AND  CORRELATION  COEFFICIENTS 

3970  DO  118  K«l,2 

3988  YAVG (K)»St /FLOAT (N) 

3990  MSO ( K ) =0 

4888  ST l K ) =0 

4818  1  10  S8(M*0 

4020 

4030  DO  120  k-1  ,2 

4848  DO  138  1*1,  N 

4858  ST (K)=ST (K) ♦ (YACT ( I ,K> -YEST ( I ,K>  >**2 

4860  138  SB ( K ) *SB(K) ♦  < YACT ( I ,K> -YAVS(K) »**2 
4878  I F  (  (1.0-ST'n/SB(Ki)  .LT.  0'  CORR(H=0. 

4888  I F  (  ( l  .0-ST  <k  )  / S B  (B  )  .•  .LT.  0)  GO  TO  125 

4090  CORR(K»-S'jRT'  1.0-ST  '  K  >  /  SB  ( F  )  > 

4188  125  I F  <  N  .EQ.  2)  bO  TO  120 

4110  STE  *K) — SORT t S T i r ) / ( N - 2 > ) 

4120  128  CONTINUE 

4130 

4140C  CALCULATE  DATA  FCc:  RE 1  URN  PEH1U0  TABLES 

4  150  DO  5  7  J=  1  ,5 

4160  F ROB* 1 . 0  - i . J/ ! LAID  Dm* RET  I  J  )  ) 

4170  IFIPROB  .LE.  0)  PROB* . 000000 1 

4100  CHS (J,l)--ALOGl-ALOblPROB)>«PHl+EPSI 

4198  CHS< J ,2) = i-ALOGl 1 .0-PROB) )**( 1.0/C) *SIGMA 

4201C  CHS(J,3)*SI5MA7/((-ALOG(PROB))**(1.0/U>) 


4210  57  CONTINUE 

4220 

4230C  CALCULATE  HEAN  SQUARE  DEVIATION  FOR  EACH  DISTRIBUTION 
4240  DO  58  I-l  ,N 

4250  21*YACT(1,1> 

4260  Z2*EPSI- (AL06 (-ALOG ( Z 1 > ) ) *PHI 

4270  Z3-BETA (2) «  C (-AL08 ( 1-Z 1 > > ** ( 1 . 0/ALPHA (2) ) ) 

42B0C  Z4*BETA(3)/( <-AL06<Zl ))**<!. 0/ ALPHA (3) ) ) 

4290  HSD<l>-HSD(l)*<22-HS<I>)*#2 

4300  HSD (2) *HSD  <2)  +  (Z3-HS(I) )  **2 

4310C  HSD(3)«HSD(3)  +  U4-HS(I)  )*»2 

4320  58  CONTINUE 

4330  HSD(l)-HSD(n/<N#PHI*«2) 

4340  HSD (2) «HSD<2) / (N*BETA(2>**2> 

4350C  HSD (3) *HSD (3)/(N*BETA(3)**2) 

4360 

4370C  CALCULATE  HEAN  AND  VARIANCE  FOR  EACH  DISTRIBUTION 

4380  HEAN ( 1 ) -EPSI +EUL£R*PHI 

4390  VAR(1)-1.6449341*PHI*»2 

4400  PARA*1 .0+1.0/C 

4410  CALL  6AHHA ( PARA , NHE ) 

4420  HEAN (2 ) "S I BHA*WHE 

4430  FACl-SIGHAt*2»WHE*»2 

4440  PARA-1. 0*2. 0/C 

4450  CALL  6AHHA (PARA,MV2) 

4460  FAC2*SIGHA*#2*WV2 

4470  VAR (2) -FAC2-FAC 1 

4480C  PARA-1. 0-1. 0/U 

4490C  CALL  GAMMA < PARA , HPC > 

4500C  HEAN(3)-SIGHA2*HPC 

4510C  PARA-1. 1-2. t/U 

4S20C  CALL  GAHHA (PARA , HP D I 

4530C  VAR(3)-SIGHA2**2*HPD-HEAN(3)**2 

4540 

4550C  WRITE  OUT  THE  DATA  FOR  EACH  DISTRIBUTION 
4560  WRITE 1 6 , 136) 

4570  136  FORHAT ( 1H1 ) 

4580  UR  I TE ( 6 , 1 35)  T I TLE , VAR  I  ABLE 

4590  135  FORHAT (///,16X,A26,A8,///> 

4600  DO  150  K-1,2 

4610  STDEVIIO-SQRT  (VAR  iK)  ) 

4620  WRI TE ( 6 , 1 60  >  IFLAG(K) ,DEF,FORH(K) 

4630  160  FORHAT <15X,A30,//, IX, A17.2X,A34) 

4640  I F  (  f  .EQ.  1)  UR  I TE ( 6 , 1 59 )  EPS  I , PH  I 

4650  159  FORHATIIX," EPSI*“6X1F10.3,/,1X," PHI*“,7X,F10.3) 

4660  IF (  K  .GT.  1)  URITE!6,161)  ALPHA (K )  , BETA (K ) 

4670  161  FORHAT(1X,"ALPHA*",6X,F10.3,/,1X,”BETA=",6X,F10.3) 

4680  URITEI6.162)  MEAN ( K ) , VAR <K> ,STDEV (K> 

4690  162  FORHATtlX, "HEAN-", 6X,F10. 3,/, IX, "VARIANCE-", 2X,F10. 3, 

4700  6/,  IX,  "STD.  DEV.  *  ",2X,F7.3> 

4710  IF  (  LOGIC  ,E0.  N  >  GO  TO  171 

4720  DO  170  1*1, N 

4730  DUH1 < I ) -YACT  > I ,  K ) 


CIO 


f  '■ 

4748  178 

DUM2 ( I ) «YEST ( I  ,K> 

• 

4738 

L2-N 

4768 

CALL  RESIDUAL (HS,DUM1 ,  DUM2 , L2) 

4778  171 

WRITE (6, 163)  CORR(K) ,ST (K) 

4788  163 

FORMAT!/, IX, "NON-LINEAR  CORRELATION  IS* , 5X ,F1 0. 7 , / 

4798 

6  , IX, "SUM  SQUARE  RESIDUALS  I S" , 6X , FI  1 . 7) 

4888 

IF (  N  .EQ.  2  )  80  TO  167 

•4;* 

4818 

WRITE (6, 164)  STE(K) 

4828  164 

FORMAT (IX, "STANDARD  ERROR  IS” , 13X ,F10. 7) 

4838  167 

WRITE (6, 166)  HSD(K) 

4848  166 

FORMAT (IX, "MEAN  SQUARE  DEVIATION  IS" ,6X ,F10. 7, ///) 

*  t 

4838  287 

NRI TE (6 , 288)  VARIABLE 

■ '» 

4868  288 

FORMAT (7X , "RETURN  PERIOD  TABLE* ,/, 6X ," YEAR" , 13X , AB> 

4878 

DO  211  Jal  ,5 

4888 

WRITE  (6,212)  RET.(J)  ,  CHS  ( J  ,  K ) 

4898  212 

FORMAT (1X,F9.2,8X,F9.2) 

4988  211 

CONTINUE 

f 

►  V 

4918 

WRITE (6, 165) 

*  t 

4928  165 

FORMAT (////) 

4938  158 

CONTINUE 

' •  Vi 

4948 

RETURN 

* 

4958 

END 

4968 

•<> 

4978 

&S 

4988 

■*!< 

4998 

■;v! 

S888C  SUBROUTINE  TO  PUT  NUMBERS  IN  ORDER  BY  ASCENDJN6  X 

■ 1 

5818 

SUBROUTINE  ORDER ( X , N) 

5828 

DIMENSION  X (N) 

5838 

INTE6ER  X , TX 

V,‘ 

I*.4 

5848 

DO  20  K»2,N 

^4*4 

5858 

J»N-k>2 

,'»V 

5868 

DO  10  1*1, J-l 

5878 

I F (  X ( I )  .LT.  X ( I ♦  1 ) )  60  TO  10 

3880 

T  X  *  X  ( I ) 

:•/ 

5898 

X(I)*X(I+i) 

'«* 

5188 

X  C I ♦ 1 ) -TX 

;'V 

5110  18 

CONTINUE 

5128  28 

CONTINUE 

5138 

RETURN 

5148 

END 

*v 

.<■*  * 

3158 

.A 

5168 

A 

3170 

518BC  SUBROUTINE  TO  HELP  PRINT  OUT  DATA 

3190 

SUBROUTINE  RES I DUAL ( X , Y ACT , YEST , N) 

r .  M 

5288 

DIMENSION  X (N) ,Y ACT (N) , YEST (N) ,DIFF (200) 

■v 

5218 

INTE6ER  X 

5228 

SSR*0 

1  i*1 

•»  y 

5238 

DO  10  1=1, N 

' , 

5240 

DIFF ( I>» (YACT ( I > -YEST< I ) )*#2 

3258  18 

SSR=SSR*DIFF ( I ) 

,  ■  * 

■ 

5268 

WRITE  <6, 15) 

Cll 

* 

LI 

5270  15  FORMAT!//, IX,"  XVALUE  YVALUE  YEST  DIFF  ",/,) 

5280  DO  25  I-1,N 

5290  MRITE<6,20)  X ( I ) , YACT ( I ) , YEST (I ) ,SQRT (DIFF  ( I ) ) 

5300  WRITE(8,20)  X ( I ) , YACT ( I ) , YEST ( I ) , SORT ( DIFF ( I ) ) 

5310  20  FORMATUX,  111,  Fll.  4,  FI  1. 4,  FI  1. 4,/,) 

5320  25  CONTINUE 

5330  RETURN 

5340  END 

5350 
5360 

S370C  SUBROUTINE  TO  EVALUATE  THE  6AMMA  FUNCTION 
S380C  PROGRAM  ADJUSTS  ALPHA  TO  BE  BETWEEN  1.0  AND  2.0 
5390C  AND  THEN  MULTIPLIES  BY  6F  TO  COMPENSATE 
5400  SUBROUTINE  6AMMA (ALPHA , AREA) 

5410  DOUBLE  PRECISION  C ( 25 > , SUM 

5420  6F-1 .0 

5430  IF (ALPHA)  1,2,3 

5440 

5450  2  PRINT, 'TROUBLE  IN  GAMMA' 

5460  AREA-1.0 

5470  60  TO  200 

5480 

5490C  FOR  GAMMA  OF  A  POSITIVE  NUMBER 
5500  3  M-INT (ALPHA) 

5510  EPSI -ALPHA-FLOAT (M) 

5520  IF (  M  .EQ.  0)  GF=GF/ALPHA 

5530  IF (  M  .EQ.  0)  ALPHA-ALPHA+ 1 . 0 

5540  IF (  M  .EQ.  0)  GO  TO  100 

5550  IF (  M  .EQ.  1)  GF-1.0 

5560  IF'.  M  .EQ.  1)  GO  TO  100 

5570  DO  10  1-2, M 

5580  10  6F«6F*(FLQAT ( I - 1 ) +EPSI ) 

5590  ALPHA-1. 0+EPSI 

5600  60  TO  100 

5610 

5620C  FOR  GAMMA  OF  A  NEGATIVE  NUMBER 
5630  1  H-I NT (ALPHA) 

5640  EPSI-ALPHA-FLOAT (M) 

5650  DO  20  1-1,2-M 

5660  J-M+ ( I - 1 ) 

5670  20  GF-GF/ (EPS I ♦FLOAT < J > ) 

5680  ALPHA-EPSI+2. 0 

5690 

S700C  COEFFICIENTS  FOR  SERIES  EXPANSION  OF  THE  GAMMA  INTEGRAL 

5710C  SEE  HANDBOOK  OF  MATHEMATICAL  FUNCTIONS  BY  ABRAMQWITZ  AND  SEGUN 

5720  100  C(l)-1. 0000000000000000 

5730  C(2)«. 5772156649015329 

5740  C(3)«-. 6558780715202538 

5750  C(4>  —  .0420026350340952 

5760  C<5)«. 1665386113822915 

5770  C(6)— .0421977345555443 

5780  C<7)—  .009621971527887 

5790  C (8) -. 007218943246663 


5800 

C(9)» 

-.0011651675918591 

5810 

C  ( 10) 

■-.0002152416741149 

5820 

C(ll) 

-.0001280502823882 

5830 

C  <  12) 

■-.0000201348547807 

5840 

C  ( 13) 

■-.0000012504934821 

5850 

C(14> 

-.0000011330272320 

5860 

C  ( 15) 

—  .0000002056338417 

5870 

C  ( 1 6) 

■6. 1 1 6095E-09 

5880 

CU7) 

■5. 0020075E-09 

5890 

CllB) 

—  1.  1812746E-09 

5900 

C  ( 19) 

» 1 . 043427E- 1 0 

5910 

C  (20) 

■7. 7823E-12 

5920 

C  (2 1 ) 

■-3. 69680E- 12 

5930 

C  (22) 

■5. IE-13 

5940 

C  (23) 

— 2.06E-14 

5950 

C  ( 24) 

— 5.4E-15 

5960 

C  (25) 

■1.4E-15 

5970 

5980C 

SUM  SERIES 

S99i  sun>e.0 

6000  DO  50  K « 1 ,25 

6010  SUM«SUH*C(K)*(ALPHA**K) 

6020  50  AREA*BF/SUH 

6030  200  RETURN 

6040  END 

60504: EXECUTE 

60604: LIMITS: 30, 100K 

60704: FILE : 07 , X7R, 5L , NEW , STRMF ILE 

60804: FIL£:08,X8R,5L,NEW,D1STFILE 

60904: FILE: 09, X9R , 5L , NEW , AB3DST30 

6 1 004 : END J  OB 


ANALYSIS  OF  STORM  DURATION 
DATA  FILE  NA6SHEAD ,  NORTH  CAROLINA 


STORM  NO 

.  DATE/TIME  OF  PEAK 

DURATION  H>350 

PEAK  H 

PEAK  T 

PEAK  DIR 

1 

5601 1013 

27 

449 

11 

98 

2 

36092718 

6 

379 

10 

101 

3 

56102703 

3 

461 

11 

83 

4 

56102800 

9 

405 

11 

75 

5 

36103100 

3 

377 

8 

69 

6 

56103109 

6 

364 

11 

74 

7 

38102103 

3 

352 

10 

91 

B 

58102118 

9 

488 

10 

75 

9 

58102212 

27 

508 

11 

88 

10 

60020100 

9 

386 

10 

93 

11 

60020118 

3 

351 

10 

101 

12 

60103121 

3 

354 

10 

97 

13 

61102421 

6 

387 

11 

99 

14 

62030721 

18 

459 

11 

114 

15 

62030909 

30 

591 

13 

94 

16 

62112806 

84 

466 

12 

97 

17 

62120200 

3 

351 

9 

88 

IB 

62120212 

15 

371 

9 

88 

19 

63020421 

6 

363 

9 

88 

20 

64092221 

21 

391 

10 

111 

21 

66061306 

3 

403 

10 

105 

22 

68011121 

12 

400 

10 

100 

23 

68022512 

15 

385 

10 

116 

24 

69022109 

3 

355 

10 

114 

23 

69030306 

3 

360 

11 

101 

26 

70102718 

12 

364 

10 

103 

27 

72052700 

15 

365 

10 

94 

28 

73021112 

33 

465 

11 

111 

29 

73021300 

3 

372 

10 

101 

30 

73022806 

9 

379 

10 

107 

31 

73120906 

3 

352 

9 

66 

32 

73012118 

6 

389 

10 

106 

33 

75070106 

9 

399 

10 

112 

34 

75070215 

9 

430 

11 

104 

33 

75112415 

9 

375 

10 

107 

36 

75112506 

3 

397 

10 

102 

1.80  STORMS  PER  YEAR 

NO.  RECORDS  WHERE  H  >  330 

=•  149  (  0.3X 

OF  58440 

RECORDS) 

MIN. 

PEAK  H  ■  351  MAX. 

■  591  MEAN  « 

401.2 

STD.  DEV. 

=  53.4 

MIN. 

DURATION  *  3  MAX 

.  *  84  MEAN  » 

12.2 

STD.  DEV 

.  =  14.9 

THE 

DATE/TIME  IS  YRMODYHR 

,  DURATION  IS  IN 

1  HOURS, 

H  ^HEIGHT)  IS  IN  CM 

T  (PERIOD)  IS  IN  SEC  AND 

DIRECTION  IS  IN 

DE6REES 

RELATIVE 

TO  THE  SHOREL 

Cl4 


I 

I 

i 


I  LEAST  SQUARES  RESULTS  ~  DURATION 

t 


EXTREMAL  TYPE  I 


F ( X ) *PR ( X<  X ) * 
EPS  I » 

PHI  * 

MEAN* 

VARIANCE* 

STD.  DEV.  = 


EXP<-EXP(-(X-£PSI)/PHI) ) 
3.918 
15.246 
12.718 
382.333 
19.553 


NON-LINEAR  CORRELATION  IS  0.8720683 
SUM  SQUARE  RESIDUALS  IS  0.6796928 
STANDARD  ERROR  IS  0.1413894 
MEAN  SQUARE  DEVIATION  IS  0.3527218 


RETURN  PERIOD  TABLE 


YEAR 

5.00 

10.00 

25.00 

50.00 

100.00 


DURATION 

36.53 

47.55 

61.78 

72.44 

83.05 


\ 

i 

> 


WE  I  BULL 


FIX) =PR ( X<  X  )  * 
ALPHA* 

BETA* 

MEAN* 

VARIANCE* 

STD.  DEV.  * 


1-EXP  I- (X/ BETA) t# ALPHA) 
1.156 
12.636  , 

12.007 

108.437 

10.413 


;  NON-LINEAR  CORRELATION  IS 

j  SUM  SQUARE  RESIDUALS  IS 

STANDARD  ERROR  IS 
MEAN  SQUARE  DEVIATION  IS 


0.9607069 

0.2166227 

0.0801878 

0.3858081 


RETURN  PERIOD 

TABLE 

YEAR 

DURAT 

5.00 

24.96 

10.00 

31.64 

25.00 

40.  15 

50.00 

46.40 

100.00 

52.52 

C15 


LEAST  SQUARES  RESULTS  -  PEAK  H 


EXTREMAL  TYPE  I 


F ( X ) -PR ( X<  X ) 
EPS  I  * 

PHI- 
MEAN- 
VARIANCE* 
STD.  DEV.  * 


i  EXP(-EXP(-(X-EPSI>/PHD) 

374.975 
48.460 
402.947 
3862.962 
62.153 


NON-LINEAR  CORRELATION  IS 
SUM  SQUARE  RESIDUALS  IS 
STANDARD  ERROR  IS 
MEAN  SQUARE  DEVIATION  IS 


0.9629306 

0.2064946 

0.0779318 

0.0998112 


RETURN  PERIOD  TABLE 


YEAR 

5.30 

10.00 

25.00 

50.00 

100.00 


PEAK  H 
478.63 
513.66 
558.90 
592.77 
626.49 


ME1BULL 


F  (  X ) -PR  (  X  <  X )  = 
ALPHA* 

BETA* 

MEAN- 

VARIANCE* 

STD.  DEV.  » 


1-EXP(-(X/DETA) ** ALPHA) 
7.888 
426. 3B8 
401.273 
3638.859 
60.323 


NON-LINEAR  CORRELATION  IS 
SUM  SQUARE  RESIDUALS  IS 
STANDARD  ERROR  IS 
MEAN  SQUARE  DEVIATION  IS 


0.9038882 
0.5192849 
0.  1235843 
0.0047993 


RETURN  PERIOD  TABLE 


YEAR 

5.00 

10.00 

25.00 

50.00 

100.00 


PEAK  H 
4.71.  13 
487.80 
505. 13 
515.95 
525.41 
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APPENDIX  D 


SPSS  COMMAND  FILE  AS  APPLIED  IN 
THE  REGRESSION  ANALYSIS 


Command  File  Listing.... 
Sample  Output  (Excerpts) 


COMMAND  FILE  LISTING  -  SPSS  REGRESSION  ANALYSIS 


1099S,T,J 

20# : I  DENT i R0CDDPS,OPSMI  TH 

309: SELECT: SPSS /SPSS 

409: SYSOUT: 43 | NULL 

509: LIMITS:  ,60K 

609: INCODE: IBMF 

70RUN  NAME: DURATION  ANALYSIS 

80VARIABLE  L 1ST: DUR , H , T , D 

90INPUT  MEDIUM: DISK 

100INPUT  FORMAT: FIXED (33X ,F4. 0, 12X,F3. 0,6* ,F2.0,6X ,F 3.0) 

I 10N  OF  CASES: UNKNOWN 

120VAR  LABELS: DUR  DURATION/H  PEAK  H/T  PEAK  T/D  PEAK  r i R / 

1 30COHPUTE: HSQ»Ht*2 
1 40COMPUTE : TSQ*T**2 
150COMPUTE: STP*H/ (981*TSQ) 

160COMPUTE: SEV»156. 13*HSQ»TSQ 

170VAR  LABELS: HSQ  H##2/TSQ  T**2/STP  H  OVER  gT**2/SEV  LH**2 
I 80REBRESSION: VARIABLES" DUR ,H,HSQ,T,TSQ/ 

190: : REGRESS  I QN*DUR  WITH  H , HSQ, T , TSQ < 1 )  RESID>0/ 
210STATISTICS: ALL 
220READ  INPUT  DATA 
230REGRESSION;VARIABLES«DUR,SEV/ 

240: : REGRESSION* DUR  WITH  SEV(l)  RES ID»0/ 

2 60S TAT I  ST  ICS: ALL 
270RE8RESSION: VARI ABLES-DUR , D/ 

280: : REGRESS  I QN»DUR  WITH  D(l)  RES  I D-0/ 

300STATIST ICS: ALL 

31 0RE6RESSI ON: VARIABLES" DUR , S  T  P / 

320: :RE6RESSI0N*DUR  WITH  STP(l)  RES  I D-0 / 

340STATISTICS:ALL 

3S0SCATTERGRAM: DUR  WITH  H , T , HSQ , STP , SEV 

360STATISTICS: ALL 

370FINISH 

3809: DATA: 08 

39099SELECT (P05DST50) 

4009 : END J  OB 


VALUE  Of  eo.oojno  IS  PR  I  NT  f  0 
A  COEFFICIENT  CANNOT  OF  COMPUTED 


00**1/  0**011  00*011  Ol'fOl  01*11*  OO’OO*  0**00*  00*11*  01*0**  0**10*  00*00* 


a7425?t.M91784S2.0O**14})t.0000SO?40.0n04**14  4.M0*2  2  04t.001]S7*S2.OO17*JtS*.nO222*7t0.M2«tS6*4.M]101S**.0n 


